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ABSTRACT 

This  paper  reports  an  investigation  of  the  dynamic  behavior 
of  a  pin-ended  inelastic  column  subjected  to  a  half-sine  dis- 
placement pulse  at  one  end  applied  with  small  eccentricity,, 
The  rates  of  end  displacement  are  such  that  the  effects  of 
axial  inertia  are  considered  negligible,,   The  column  is  re- 
placed by  a  lumped  parameter  mathematical  model  and  the  equa- 
tions for  the  model  are  solved  by  a  high  speed  digital  com- 
puter.  The  stress-strain  relation  during  loading  is  assumed 
to  be  of  the  Ramberg-Osgood  type„   Unloading  occurs  along  a 
line  with  slope  equal  to  that  of  the  tangent  at  the  origin,, 
Slendemess  ratios  of  50,  70,  and  100  and  two  different  values 
of  eccentricity  are  considered  in  the  investigation.  The  max- 
imum load  supported  by  the  column  and  the  residual  lateral 
deflection  at  the  center  of  the  column  are  found  as  functions 
of  the  displacement  pulse  height  and  of  the  displacement  pulse 
duration. 

It  Is  shown  that  an  inelastic  column,  like  the  elastic 
column,  can  support  dynamic  loads  considerably  in  excess  of 
the  static  load  capacity.  It  Is  also  shown  that  within  the 
range  of  pulse  durations  studied,  the  smallest  residual  de- 
flections generally  accompany  the  largest  maximum  loads. 
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1.    INTRODUCTION 

In  the  study  of  structures  the  column  ranks  as  one  of 
the  more  provoking  elements  as  judged  from  the  amount  of 
literature  published  on  Its  behavior.   The  development  of 
elastic  column  theory  is  comparatively  advanced  and  dates 
back  to  1757  when  Suler  published  his  analysis  of  the  strength 
of  columns.   In  more  recent  years  inelastic  column  theory  has 
been  developed  to  the  point  where  a  satisfactory  explanation 
has  been  obtained  for  the  phenomenon 0   Shanley  /l/*  has  been 
credited  with  revising  and  reinterpreting  inelastic  column 
theory  from  previously  accepted  erroneous  conjectures.   Util- 
izing studies  of  a  highly  idealized  column  he  concluded  that 
the  tangent  modulus  load  is  the  maximum  load  at  which  an  in- 
itially straight,  centrally  loaded  column  will  remain  straight; 
loading  beyond  this  load  causes  bending  which  produces  per- 
manent deformation. 

The  dynamic  behavior  of  structures  has  become  an  increas- 
ingly important  engineering  consideration  in  recent  years  and, 
with  the  aid  of  more  sophisticated  devices  for  the  problem 
solution,  many  more  phenomena  can  be  Included  In  an  analysis 
than  could  previously  be  reasonably  considered.   The  behavior 
of  dynamically  loaded  elastic  columns  has  received  consider- 
able attention  and  a  number  of  investigations  have  resulted. 

Hoff  and  associates  /2, 3/  considered  the  case  of  an  elas- 
tic column  of  uniform  cross  section,  initially  curved  in  the 

*  Numbers  in  /  /  refer  to  entries  in  the  bibliography 


shape  of  a  half-sine  wave,  subjected  to  constant  velocity 
loading  such  as  that  encountered  in  compression  tests  in 
commercial  testing  machines.   It  was  shown  that  the  maximum 
loads  supported  by  the  column  are  greater  than  the  Suler  load 
when  the  loading  is  rapid,  the  column  is  slender  and  the  ini- 
tial deflections  are  small.   Sevin  /4/  confirmed  the  previ- 
ous results  while  including  the  effects  of  axial  inertia 
which  were  not  considered  by  Hoff.   He  Inferred  that,  for  an 
elastic  column,  axial  inertia  effects  are  of  negligible  im- 
portance with  regard  to  the  gross  behavior  of  conventional 
structural  columns,  regardless  of  the  initial  deflected  shape, 
end  constraint,  or  type  of  axial  loading. 

Taylor  /5/  investigated  the  case  of  an  eccentrically 
loaded  elastic  column  subjected  to  a  half-sine  force  pulse 
and  he  Included  the  effects  of  axial  inertia.   Using  a  limit- 
ing value  of  extreme  fiber  strain  as  the  failure  criterion 
he  determined  the  column  failure  load  quantitatively.   He 
demonstrated  that  elastic  columns  will  support  loads  much 
greater  than  the  Euler  load  when  subjected  to  rapid  loading 
and  that  axial  Inertia  effects  become  significant  with  this 
type  of  loading  as  the  rapidity  of  loading  is  increased. 

Gerard  and  Becker  /6/  studied  column  behavior  under  con- 
ditions of  impact,  by  utilizing  the  properties  of  a  propa- 
gated elastic  compressive  stress  wave  produced  on  failure 
of  a  tension  specimen.   They  concluded  that  a  column  can 
momentarily  support  a  dynamic  compressive  stress  of  any 


magnitude  and  that  buckling  may  be  confined  to  a  length  less 
than  the  entire  length  of  the  column , 

The  inelastic  column  was  investigated  by   Chav;la  /7/ 
in  essentially  the  same  manner  as  Hoff  (with  an  initial 
curvature  or  effective  eccentricity)  and  he  found  that  for 
relatively  slow  rates  of  end  displacement  the  maximum  column 
load  was  less  than  the  tangent  modulus  load.   Brooks  and 
Wilder  /&/   studied  the  effects  of  dynamic  loading  on  an 
Idealized  H-section  Inelastic  column  with  an  initial  lateral 
deflection  using  a  constant  velocity  end  displacement.   From 
their  results  they  concluded  that  the  static  maximum  load 
may  be  employed  as  a  conservative  estimate  of  the  maximum 
load  of  a  column  regardless  of  the  rate  of  end  displacement 0 

Investigations  of  the  dynamic  behavior  of  inelastic 
columns  are  limited  in  number  and  the  problem  of  determin- 
ing the  behavior  of  inelastic  columns  subjected  to  pulse 
loading  appears  to  have  been  neglected,,   Whereas  constant 
velocity  loading  leads  to  the  ultimate  failure  of  the  column, 
pulse  loading  may  not  necessarily  result  in  failure,,   Pulse 
loading  allows  the  column  to  recover  from  a  deflected  posi- 
tion with  an  amount  of  residual  deflection  while  retaining 
a  capacity  to  sustain  subsequent  load.   The  problem  con- 
sidered in  this  investigation  is  that  of  determining  the 
dynamic  behavior  of  a  pin-ended  inelastic  column  subjected 
to  a  half-sine  displacement  pulse  at  one  end  applied  with 
small  eccentricity.   The  problem  involves  the  prediction 


of  the  maximum  force  developed  In  the  column  for  a  specified 
pulse  height  and  duration  of  loading.   The  problem  also  in- 
cludes the  prediction  of  the  residual  lateral  deflection 
resulting  from  the  pulse  loading.   The  effects  of  pulse  height 
and  duration  are  investigated  for  tv/o  values  of  eccentricity 
ana  three  column  slenderness  ratios 0 


2 •    DEVELOPMENT  OF_  THE  EQUATIONS  FOR  THE  MATHEMATICAL 

MODEL 

The  column  considered  in  the  investigation  is  initially 
straight,  pin-ended,  and  of  constant  cross-section,,   The  ef- 
fects of  axial  inertia  and  shear  strains  are  neglected,  but 
non-linear  axial  strain  components  due  to  bending  at  the 
center  portion  of  the  column  length  are  considered 0   The 
effect  of  a  non-linear  stress-strain  curve  on  loading  and 
departure  from  this  curve  on  unloading  is  also  included,, 

The  loading  on  the  column  is  produced  by  axial  displace- 
ment of  the  pin  at  one  end.   The  loaded  end  is  free  to  rotate 
about  the  pin,  which  is  eccentrically  placed.   The  other  end 
is  constrained  by  a  fixed  pin  of  equal  eccentricity c   The 
end  displacement  is  Imposed  by  a  moving  ram  executing  a  half- 
sine  displacement  pulse0   The  ram  can  exert  only  compressive 
force  on  the  column  (see  Fig  0  l)o 

The  problem  is  developed  and  solved,  not  in  terms  of  the 
real  column,  but  in  terms  of  the  lumped  parameter  mathematical 
model  shown  in  Flg0  1„   The  equations  for  the  model  are  then 
solved  in  a  high  speed  digital  computer,,   The  force~def lection 
characteristic  of  the  springs  is  derived  from  the  material 
stress-strain  diagram,,   The  assumed  behavior  is  Indicated  in 
Fig.  2  where  the  solid  curve  OABD  represents  the  variation 
of  compressive  stress  O    under  montonically  increasing  com- 
pressive strain  £  „   If,  upon  reaching  point  B,  the  strain- 
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Fig.   1,   Geometry  of  the  Mathematical  Column  Model. 


Fig.   2.   Stress  Behavior  In  the  Column    Model  Springs, 


ir.g  is  reversed,  the  path  will  be  along  BC  whose  slope  is 
that  of  the  tangent  to  the  solid  curve  at  G.   A  second  re- 
versal will  cause  the  compressive  stress  to  increase  again 
following  CB0   After  again  reaching  the  curve  at  B,  the  path 
will  be  along  the  solid  curve  toward  D0   Additional  strain 
reversals  would  result  in  additional  "detours"  parallel 
to  BC0   In  the  study  the  stress-strain  relationship  after 
departure  from  the  compressive  stress- strain  curve  is  con- 
sidered linear  into  the  tensile  stress  region.   The  maximum 
tensile  stresses  in  the  column  model  are  limited  to  values 
near  the  nominal  yield  stress  of  the  material „   It  is  pre- 
sumed that  any  tensile  stress  in  the  column  model  springs 
is  insufficient  to  cause  plastic  straining0 

This  behavior  is  incorporated  in  the  springs  by  taking 
the  stiffness  of  each  as  K=A,  ET/JL ,  where  A,  is  half  the  col- 
umn cross-sectional  area,  ET  is  the  tangent  modulus,  and  Si 
is  the  unstrained  length  of  a  springo   Again  referring  to 
Fig.  2,  at  any  point  on  the  solid  curve  Er  is  the  slope  of 
the  curve.   For  points  below  the  curve  E~=E  (Young's  Modulus), 

This  model  having  only  two  degrees  of  freedom  was  cnosen 
for  the  investigation  in  order  to  limit  the  computer  solution 
time  to  reasonable  values,  but  yet  to  retain  the  effects  of 
Inelastic  column  behavior,   It  is  felt  that  this  model  is  a 
good  compromise  in  that  it  approaches  the  deflected  shape 
of  the  inelastic  column  and  also  retains  the  effects  of 
lateral  inertia  and  of  the  longitudinal  strains  on  each  side 


of  the  column  due  to  both  the  axial  load  and  bending. 

From  the  free  body  diagram  of  ha]f  the  column  model 

shown  in  Fig.  3,  using  d,Alembert,a  principle,  the  follow, 

lng  equations  may  be  written: 


?  =  f;  +  rz 


(i) 


Q=^- 


2-  dtz 


(2) 


P  (e  cos-e-+  -^-  sm-e-)-  Q(-^-  cos-e-  -  e  sine)  = 

(F.-Fj^cose 

Dividing  the  latter  by  cos-0-  and  substituting  for  Q  from 
equation  (2) 


P(e  +  -^-tan-e-)  - 


(3) 


Fig.  3o  Free  Body  Diagram  of  Half  the  Column  Model. 
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From  kinematic  considerations  of  the  column  model  in 
a  deflected  position  the  following  equations  may  be  written: 

cf,  =  <£m   +  2L^sin-e  (4) 

6Z=  eta  -a^sin-e-  (5) 

where  o,  and  ^  are  the  axial  shortening  of  Spring  1  and 
Spring  2  respectively  and  Om  is  the  average  axial  short- 
ening at  the  center  of  the  model.  Further 

x  =  eta  +  2.*4f-(j-cos-e-)  +  z  e  sin-e-       (6) 

y   =T"  5ine  ""  e(l-COS-e)  (7) 

where  x  is  the  end  displacement  at  the  loaded  end  of  the 
model  and  y  is  the  lateral  deflection  at  midlength  in  the 
model. 

The  following  approximations  are  consequences  of  assum- 
ing that  -©*  is  small  and  e/L  is  small: 

~  »  e  t  a  n  -©■ 

—  sm-e  »  e0~cos_^) 

t  an -0  «  sin-e-  ^  O 

i  -  cose  ^  "k"^"2. 

An  evaluation  of  the  errors  resulting  from  these  approxima- 


tions  is  given  in  App.  VI. 

Equation  (7)  now  becomes: 


y=^^ 


(8) 


Now,  substituting  the  approximations,  together  with  equa- 
tion (8),  into  equations  (3)  through  (6)  and  rearranging, 
results  in 


Solving  for  dm    and  substituting  in  the  above  expressions 
for  o,  and  cSg.  gives 

4  =  *  - -ty*  -  *£ y  +  4§y  (10) 

4-*- •£/-¥?-¥*  (id 

Since  the  model  has  but  two  degrees  of  freedom,  only  a 
limited  correspondence  may  be  established  between  column  and 
model  parameters .   Energy  considerations  are  used  for  this 
purpose. 
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To  determine  this  correspondence  both  the  column  and  the 
model  are  considered  for  the  elastic  case.   Figure  4  indicates 
an  elastic  column  in  a  deflected  position  during  free  vibra- 
tion.  It  has  an  assumed  half- sine  deflected  shape  and  the 
motion  is  considered  sinusoidal  with  time. 

The  maximum  potential  energy  is  found  from  the  elastic 
strain  energy  due  to  bending  and  is 

.2- 


Potentlal  Energy  =  -r- 


^MW^ 


v.'here  4- 4   is  "the  curvature  of  the  axis  of  the  column  and  I 
is  the  moment  of  inertia  of  the  column  cross-sectional  area, 
I  is  equal  to  o^A,  where  A  is  the  column  cross- sectional 
area.  This  can  be  shown  to  be 

Potential  Energy  =  """*  ^  *  jf^  (column) 

The  maximum  kinetic  energy  of  the  column,  neglecting 
axial  inertia  effects,  is 


Kinetic  Energy  =-|  )AU^;  \  yad*. 


which  can  be  shown  to  be 

Kinetic  Energy  =4  ^  ^n  V1^  (column) 

where  m.  is  the  mass  per  unit  length,  Gon  is  the  circular 
frequency  of  first  mode  lateral  vibration  of  the  column  and 
L  is  the  column  length.   Setting  the  maximum  potential  energy 
equal  to  the  maximum  kinetic  energy  and  solving  for  the  circular 
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y  -  y*  sin  ~ 


Fig.  4.  Elastic  Column  in  Free  Vibration. 


Fig.  5*  Elastic  Column  Model  in  Free  Vibration. 
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frequency  of  first  mode  lateral  vibration  of  the  column  re. 
suits  in 


^=TrV^Uf  (column) 

Figure  5  shows  the  model  with  non-eccentric  pin  connec- 
tions in  a  deflected  position.  The  stiffness  of  each  spring 
iS  oi*t  and  the  motion  of  the  mass  is  again  considered  sinusoi- 
dal with  time.  The  maximum  potential  energy  as  found  from 
the  elastic  strain  energy  in  the  springs  is 

Potential  Energy  =  £  -M=- (  c(,  ■+  <£z   ) 

Considering  the  model  to  be  in  free  vibration  without  axial 
load  the  average  spring  deflection  is  zero  or  o,  =  -  o 2  • 
Under  these  conditions  it  can  be  shown  from  equations  (10) 
and  (11)  that 

and  the  model  maximum  potential  energy  becomes 

Potential  Energy  =  — \|t' y^  (model) 

The  maximum  kinetic  energy  of  the  model  neglecting  axial 
inertia  effects  is 

Kinetic  Energy  =  •£  m  OJ^  y^  (model) 

Now,  setting  the  maximum  potential  energy  equal  to  the  max- 
imum kinetic  energy  and  solving  for  the  circular  frequency 
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of  first  mode  lateral  vibration  of  the  model  results  in 


^  V  ^AL^  (model) 

To  determine  the  undeflected  length  £  of  the  model 
springs  the  potential  energies  of  the  column  and  model  are 
equated.   Taking  the  center  deflections  ym  of  column  and 
model  to  be  equal,  X    is  found  to  be 


TT 


4  »-  (12) 


To  determine  the  relation  between  column  and  model 

masses  the  kinetic  energies  of  the  column  and  model  are 

equated.   Again  taking  the  center  deflections  y^  to  be  equal, 
m  is  found  to  be 

*-*£■  (13) 

With  the  two  previous  correspondences  it  can  be  shown  that 
the  circular  frequencies  of  the  column  and  the  model  are  the 
same. 

Now,  considering  the  differential  equation  of  motion 
(as  previously  derived  for  the  model)  for  the  elastic  case, 
together  with  the  results  of  the  elastic  energy  development, 
Euler  load  and  end  displacement  correspondences  between  the 
column  and  model  are  determined.   These  will  be  used  later  in 
the  analysis  to  define  non-dimensional  load  and  displacement 
parameters.   If  the  acceleration  and  the  eccentricity  of  the 
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column  model  are  both  zero,  equation  (9)  becomes 

Py  =  (F,-FO^  (14) 

Equation  (14)  applies  to  static  buckling  of  an  elastic  column 
with  zero  eccentricity.   Now,  for  the  elastic  case 

(r=^&a,  ana        f,=  ^^ 

From  equations  (10)  and  (11) 
Thus 


r  _r  -Hi  pA,E  v 


Substituting  this  into  equation  (14),  the  Euler  load  PE  of 
the  column  model  is  found  to  be 

P  -  TT4  A,E 

where  L/o  is  the  slenderness  ratio  of  the  column.   In  terms 
of  the  total  cross-sectional  area  the  equation  becomes 

The  Euler  load,  in  terms  of  the  Euler  displacement,  is 

and,  substituting  for  K  and  PE  ,  the  Euler  displacement  becomes 
X,=  ^1  (16) 
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With  these  results  determined  from  elastic  behavior  of 
the  column  and  model,  the  equations  for  the  strains  and  for 
dynamic  equilibrium  of  the  inelastic  column  model  are  develop- 
ed in  a  convenient  non-dimensional  parametric  form.   From 
equation  (10)  the  strain  in  Spring  1  is 

c»   X        3  2.   L 


or 


'32.  Li-    L3^     U-  ^  '  LT 
Defining  a  non-dimensional  lateral  deflection  as  V=  y/e  and 
making  the  substitution  y=Ve,  together  with  the  substitution 
for  the  Euler  displacement  from  equation  (16),  the  strain  may 
be  expressed  as 


tr-^-[\-i^-{ftv+\v]     (17 


Applying  this  same  reasoning  for  the  strain  in  Spring  2  re- 
sults in 

4- 


6,= 


TV 


z 


Lt-i(t^"ct)V"tv] 


80-7^  x£   aVo)v   v^;  v   ^v    (18) 


From  these  equations  for  the  strains  at  the  center  por- 
tion of  the  column  model  it  can  be  recognized  that  the  first 
terms  are  the  dimensionless  end  displacements,  the  second 
terms  are  the  corrections  to  this  end  displacement  due  to 
the  angular  rotation  of  the  rigid  rods,  the  third  terms  are 
the  corrections  due  to  the  change  of  the  effective  eccentri- 
city with  lateral  deflection,  and  the  fourth  terms  result 
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from  bending. 

The  relation  between  the  strains  and  stresses  during 
Initial  loading  of  the  column  model  is  given  by  the  Ramberg- 
Osgood  /9/  equation: 

where  £  is  the  compressive  strain,  0"  is  the  compressive 
stress,  0"u  is  the  yield  stress  of  the  material  and  n  is 
the  shape  parameter  of  the  stress-strain  curve.   By  defining 


be  written  as 


a  non-dimensional  stress  S  as  -i=e-  the  previous  equation  may 


t-s  +  ^±rsn 


(19) 


This  is  the  equation  which  is  used  to  determine  the  stress 
in  each  spring  as  the  model  is  initially  loaded.  For  the  com- 
pressive stress-strain  relation  considered  in  this  investi- 
gation the  shape  parameter  is  n=10  and  the  dimenslonless 
yield  stress  is  Su=0.004.   The  stress-strain  curve  is  shown 
in  Fig.  6  and  is  typical  of  24S-T  aluminum  alloy  material. 

As  each  of  the  springs  unloads,  it  leaves  the  stress- 
strain  curve  following  a  straight  line  of  slope  9-^- -  ]_  „ 
As  subsequent  loading  takes  place  in  either  spring  the  same 
straight  line  is  retraced  and,  as  the  stress  reaches  a  level 
that  puts  it  on  the  stress-strain  curve  again,  it  follows 
the  curve  until  further  unloading  occurs. 

The  differential  equation  of  motion  is  developed  in 
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nondlmensional  form  by  dividing  equation  (9)  by  the  eccentri- 
city e  and  the  square  of  the  circular  frequency  of  first  mode 
lateral  vibration  COn  and  substituting  from  equation  (1)  for 
P.   After-  rearranging  terms  the  equation  becomes 

I   d*C//e)  _  4 


*M  (W(*  +  'H^)£l 


^Z    d  t2-      TT' 

A  dimensionless  time  T  is  defined  as  T  =  LOnt,  so  that 
U>*  (dtf  =  (dT)1.   Then'   il^M.  becomes   ^i^/e)   . 
By  making  the  substitution  V  for  y/e  this  becomes  ——?_  which 

is  symbolized  as  V.   By  substituting  V  and  V  and  recognizing 

vr 

the  form  /  =  as  the  previously  defined  non-dimensional  stress, 

the  differential  equation  of  motion  becomes 

V  =  4,  (L/?)z  [(S.+S0CV+  0-(S.-St)#-]      (20) 

The  two  terms  in  the  brackets  can.  be  recognized,  respective- 
ly, as  proportional  to  the  external  bending  moment  due  to  the 
force  on  the  end  of  the  column,  and  the  Internal  bending  mom- 
ent due  to  the  spring  forces. 

In  order  to  solve  this  equation  it  is  necessary  to  specify 
the  forcing  function.   This  is  accomplished  by  specifying  the 
displacement  of  the  ram  (as  shown  in  Fig.  1 )  as  a  function  of 
time.   When  the  ram  is  in  contact  vlth  the  column  model  the 
ram  displacement  is  the  end  displacement  of  the  column.   If 
the  ram  is  not  in  contact  with  the  end  of  the  column,  the 
force  is  zero  (S|=-S^)  and  the  column  exhibits  free  vibration. 
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In  either  case  the  conditions  of  compatibility  and  dynamic 
equilibrium  are  maintained. 

The  ram  displacement  is  a  half- sine  pulse  of  height  X 
and  duration  T^  .   To  enable  a  selection  of  dimensionless 
parameters  having  a  relationship  to  the  column  the  follow- 
ing dimensionless  quantities  are  defined: 

(a)  the  height  of  the  half-sine  ram  displacement  pulse 
as  a  fraction  of  the  Euler  end  displacement; 

(b)  the  ratio  of  the  natural  period  Tn  of  first  mode 
lateral  vibration  of  the  column  to  the  duration  T$ 
of  the  half-sine  pulse; 

With  these  dimensionless  quantities  and  the  previously  defined 
dimensionless  time.   The  dimensionless  ram  displacement 
becomes 

X=«xsin^T        for  o  ^  T  ^  *£- 

and  (21) 

"X=  O  for  T  >  Sr 

V 

To  determine  the  maximum  force  and  residual  lateral 

deflection  resulting  from  a  specified  ram  displacement  it  is 

necessary  to  compute  the  dynamic  response  of  the  model  to  the 

loading.   This  Is  accomplished  through  a  numerical  integration 

process  applied  to  the  differential  equation  of  motion  and  a 

numerical  iteration  scheme  to  satisfy  the  relations  among  the 

strains,  stresses,  end  displacement  and  lateral  deflection. 
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3.    THE  NUMERICAL  SOLUTION 

The  equations  for  the  model  are  solved  utilizing  FOR- 
TRAN programming  and  a  Control  Data  Corporation  1604  high 
speed  digital  computer.   The  investigation  is  accomplished 
with  two  basic  programs,  COL  1  and  COL  2.   The  former  is  em- 
ployed primarily  to  obtain  graphical  output  directly  to  study 
the  dynamic  response  of  the  model  to  a  specified  loading.   Pro- 
gram COL  2  is  employed  to  determine  the  numerical  values  of 
maximum  force  and  residual  lateral  deflection  for  a  specified 
loading  on  the  model.   The  complete  programs  in  FORTRAN  lan- 
guage, together  with  the  program  notation,  are  given  in  Appen- 
dix IV. 

The  numerical  integration  procedure  is  the  same  for 
both  programs,  but  the  program,  COL  2,  employed  to  determine 
maximum  values  is  more  comprehensive.   For  this  reason  the  de- 
tailed discussion  is  limited  to  this  program.   A  discussion 
of  the  differences  between  these  programs  follows  the  detail- 
ed discussion.   Figure  7  is  a  simplified  block  diagram  of 
Program  COL  2  and  Fig.  8  is  a  simplified  block  diagram  of  the 
subroutine  used  in  both  programs  to  determine  the  stress  from 
the  strain. 

Referring  to  Fig.  7,  the  input  to  the  program  Includes 
the  following  dimensionless  parameters: 

(a)  L/(p  (slenderness  ratio)  and  e/o  (eccentricity  ratio) 
to  prescribe  the  column  model, 

(b)  ex.  (ram  pulse  height)  and  Q-  (ratio  of  natural  period 
of  first  mode  lateral  vibration  to  ram  pulse  duration) 
to  prescribe  the  forcing  function  on  the  model. 
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Fig.  7.  Simplified  Block  Diagram  of  Computer  Program  COL  2 
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(c)  n  (Ramb erg- Osgood  shape  parameter)  and  3y   (yield 
stress  of  the  material)  to  prescribe  the  compressive 
3tress-strain  curve, 

(d)  AT  (time  increment  )  to  prescribe  the  incremental 
values  in  the  numerical  integration  process. 

The  initial  conditions  are  then  set  to  start  the  numerical 

integration  process. 

Henceforth,  the  dimensional  quantities  and  the  dimension- 
less  quantities  are  referred  to  synonomously  except  where  con- 
fusion may  result,  i.e.,  lateral  deflection  for  dimensionless 
lateral  deflection,  etc. 

The  numerical  integration  scheme  used  in  the  programs 
is  a  Runge-Kutta  process  /10/.  This  process  predicts  the 
lateral  deflection  V  and  the  lateral  velocity  V  at  a  time 
T+&T  from  four  estimates  of  the  lateral  acceleration  V  cal- 
culated in  the  time  increment  AT.   Accuracy  is  of  the  order 
of  (  A  T)3  .  Within  this  predictor  loop  further  calculations 
are  made  to  satisfy  the  relations  among  the  strains,  stresses, 
end  displacement  and  lateral  deflection. 

The  ram  displacement  is  computed  from  equation  (21). 
The  values  of  ram  displacement  and  of  lateral  deflection  are 
then  used  to  determine  the  strain  in  Spring  1  from  equation 
(17).  The  main  program  then  calls  Subroutine  STRESS   Because 
this  routine  is  written  with  dummy  variable  inputs  to  accommo- 
date the  strains  and  stresses  in  both  Spring  1  and  Spring  2, 
the  following  discussion  is  presented  based  on  a  general 
strain  €  and  a  general  dimensionless  stress  S. 
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Subroutine  STRESS,  shown  in  Fig.  8,  is  entered  with  the 
present  value  of  strain  £  and  the  previous  values  of  stress  Sa 
and  strain  €a  calculated  for  the  preceding  time  T-AT.   The 
solid  curve  in  Fig.  9  represents  the  diraensionless  compres- 
sive stress-strain  curve.   The  dashed  vertical  lines  indicat- 
ed by  A,  B,  and  C  in  the  figure  are  representative  of  con- 
ditions which  must  be  considered  in  arriving  at  a  value  for 
the  stress  3  from  a  value  of  the  strain  €. 

Referring  now  to  Figs.  8  and  9,  if  the  present  value  of 
strain  is  negative  this  condition  is  represented  by  A.   The 
spring  has  unloaded  from  the  stress-strain  curve  following 
a  straight  line  with  a  slope  equal  to  unity  and  the  stress 
is  calculated  from 

S»Sa  +  U-*»)  (22) 

This  equation  is  also  used  to  determine  the  stress  S  if  €=C. 
If  £  is  positive  there  are  three  possibilities  to  be  consider- 
ed,,  These  are: 

(a)  loading  is  occurring  along  the  compressive  stress- 
strain  curve; 

(b)  unloading  is  occurring  on  a  unity  slope  line; 

(c)  loading  is  occurring  on  a  unity  slope  line. 

To  determine  which  of  these  possibilities  exists  the  incre- 
mental change  in  strain,  A£=  £-£a   ,  is  calculated.   If 
this  quantity  is  negative  this  condition  is  represented  by 
B,   The  spring  is  unloading  on  a  unity  slope  line  and  the 
stress  is  calculated  from  equation  (22).   If  the  incremental 
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Fig*  8.  Simplified  Block  Diagram  of  Computer  Program  for 
Subroutine  STRESS ♦ 
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change  in  strain  is  positive  two  possibilities  exist 0  Either 
the  stress  is  on  the  compressive  stress- strain  curve  or  the 
stress  is  Increasing  on  a  unity  slope  line0   This  situation 
is  represented  by  C  in  Fig0  9„   The  stress  level  is  deter- 
mined by  calculating  a  value  of  stress  from  equation  (22) 
and  by  calculating  a  value  of  stress  from  the  Ramb erg- Osgood 
equation  (19)  and  then  selecting  the  smaller  of  these  two 
values  as  the  stress  existing  in  the  sprlng0   If  the  incre- 
mental change  in  strain  is  zero,  the  stress  S  is  equal  to  S3 

To  calculate  the  value  of  stress  from  the  Ramberg-Osgood 
equation  it  is  necessary  to  compute  the  real,  positive  root 
of  the  equation.   This  is  accomplished  by  a  Newton-Raphson 
Iteration  /10/  where  the  strain  f  is  used  as  the  initial 
estimate  for  the  root  S*.  »     The  stress  S  could  have  been  de- 
termined from  the  strain  €  each  time  by  calculating  a  value 
from  equation  (22)  and  a  value  from  the  Ramberg-Osgood  equa- 
tion (19)  and  selecting  the  smaller  value*   However,  this 
would  involve  an  iteration  for  each  stress  calculation,  thus 
increasing  the  computer  tlme0   Therefore,  this  method  was  not 
used0 

After  S  is  determined  the  value  for  the  strain  €,  is 
calculated  from  equation  (18) .   The  main  program  again  calls 
Subroutine  STRESS  to  determine  the  value  of  S|  , 

At  this  point  the  program  calculates  the  sura  S.+Sa,, 
which  is  directly  proportional  to  the  force  on  the  end  of 
the  model,  to  determine  if  it  is  compressive  (positive)  or 
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tensile  (negative).   As  previously  discussed  for  the  model 
configuration  (Fig.  1),  the  ram  cannot  transmit  a  tensile 
force  to  the  model.   Therefore,  if  this  sum  Is  negative  the 
program  initiates  an  Iterative  interpolation  scheme  to  deter- 
mine the  end  displacement  of  the  model  which  will  cause  the 
sum  Si+Sg.  to  be  zero.   This  involves  recalculating  the   strains 
and  then  the  stresses  for  each  end  displacement  until  the  force 
is  zero.   The  model  is  then  in  a  state  of  free  vibration.   This 
free  vibration  continues  until  the  ram  again  contacts  the  column 
or  until  the  calculation  terminates. 

Once  these  states  of  stress  (end  force  either  positive 
or  zero)  are  determined,  the  lateral  acceleration  V  is  cal- 
culated from  equation  (20).   After  four  passes  are  made  to 
this  point  in  the  program  to  determine  the  values  for  the 
Runge-Kutta  predictions,  as  previously  discussed,  the  lateral 
velocity  V  and  the  lateral  deflection  V  are  calculated 
for  the  end  of  the  time  interval  A  T. 

At  this  point  the  program  computes  the  force  which  can 
be  shown  to  be 


^  =  ^(l/?)Xs1  +  s,) 


and  a  continuous  check  Is  made  to  determine  the  maximum  value 
of  this  quantity  which  is  stored  for  print-out.   This  contin- 
uous check  is  also  made  to  determine  the  maximum  value  of 
compressive  stress  in  each  of  the  springs.   These  values  are 
stored  for  use  In  strain  energy  calculations  to  be  made  later, 
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A  continuous  check  Is  also  made  to  determine  the  extreme  stress 
(tensile  or  compressive)  in  each  spring  for  print-out.   The  pro- 
gram then  calculates  the  ram  work  done  on  the  model  for  the  pre- 
ceding time  increment  and  adds  this  quantity  to  the  previous 
work  done  by  the  ram. 

This  whole  process  continues  until  the  ram  displacement 
pulse  has  terminated  and  the  column  is  left  In  a  state  of 
free  vibration.   Then  a  continuous  check  is  initiated  to 
determine  the  time  at  which  S,  =Si=0o   The  lateral  deflec- 
tion at  this  point  is  the  residual  deflection  y»-  in  the 
model,  as  the  lateral  acceleration,  the  longitudinal  end 
force  and  the  internal  moment  are  all  zero,  and  the  residual 
strains  determine  the  permanent  deformation  that  has  resulted. 

After  the  residual  deflection  has  been  calculated  a 
switch  is  set  to  detour  this  portion  of  the  program.  The 
new  path  includes  a  continuous  check  to  determine  the  time 
at  which  the  lateral  velocity  is  zero  (point  of  extreme  lateral 
deflection  and  zero  kinetic  energy).  At  this  point  the  net 
energy  input  to  each  of  the  springs  is  calculated  from  the 
area  under  the  stress-strain  curve  using  the  present  values 
of  stress  and  the  previously  stored  maximum  values  of  the  com- 
pressive stress.   This  energy  is  a  combination  of  the  strain 
energy  of  plastic  deformation  and  of  the  potential  energy  due 
to  existing  elastic  strain  in  the  model.   This  energy  should 
be  equal  to  the  ram  work  done  on  the  model. 

When  the  energy  has  been  calculated  the  program  then 
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prints  all  input  quantities  and  the  results.   The  results 
include  the  values  of  Pm/Pt  ,y,-/^  ,Slrnax  ,  Szmw     and  the 
time  T  that  each  of  these  occur,  as  v.'ell  as  the  work  and 
energy.   The  times  are  included  to  complement  the  time  plots 
from  the  other  program  as  an  aid  in  visualizing  some  of  the 
dynamic  behavior  and  to  determine  the  problem  time  termina- 
tion of  the  program. 

The  program  employed  to  obtain  graphical  output  directly 
is  Program  COL  1  (the  FORTRAN  language  listing  is  given  in 
Appendix  IV) .  This  program  differs  basically  from  Program 
COL  2  described  in  the  preceding  paragraphs  in  that  it  stores 
the  variable  quantities  for  each  time  T  in  array  form  instead 
of  storing  only  the  maximum  values.  These  array  quantities 
are  then  available  for  Subroutine  GRAPH  (not  shown  in  the 
listing  but  available  on  the  FORTRAN  Library  Tape  at  the  U.S. 
Naval  Postgraduate  School,  and  also  available  from  the  Control 
Data  Corporation).   Subroutine  GRAPH  is  designed  to  accept 
a  maximum  of  900  points  for  each  variable  to  be  plotted.  Thus 
the  program  is  written  to  calculate  only  900  points  while  main- 
taining  the  necessary  accuracy  for  the  numerical  integration 
process  compatible  with  the  time  increment  AT. 

In  addition  to  the  input  data  required  in  Program  COL  2, 
axis  scaling  and  multiple  curve  data,  graph  title  data  and 
curve  label  data  are  required  for  Program  COL  1„   In  addi- 
tion to  the  graphical  output  the  program  prints  the  input 
quantities  and  the  values  of  T,  x/x£,  6,  ,  €t  ,  S.  ,  S^,  y/^ 
and  p/PEfor  every  fifteenth  point.   As  all  of  these  quantities 
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have  been  stored  as  9C0  point  arrays  any  combination  of  these 
variables  may  be  plotted  as  a  function  of  another.   Subroutine 
GRAPH  converts  these  900  point  arrays  to  input  data  on  a  mag- 
netic tape  for  a  Control  Data  Corporation  160  computer  which 
is  programmed  to  provide  output  as  input  to  a  plotter. 

An  example  of  one  set  of  curves  used  in  the  investi- 
gation is  shown  in  Appendix  I.   The  example  Includes  three 
graphs  with  two  curves  on  each  graph.  The  first  graph  is  a 
plot  of  €,  and  £z  versus  T.  The  second  graph  is  a  plot  of 
Si  and  S^  versus  T  and  the  third  graph  is  a  plot  of  y/P  and 
F/Pg;  versus  T,   Other  examples  of  graphical  output  are  shown 
in  Appendices  II  and  III.  However,  the  graphs  shown  in 
Appendix  II,  which  are  for  very  slow  loading  approaching  the 
static  case,  have  beer  obtained  with  a  modified  version  of  the 
program.  The  modification  includes  calculating  1800  points 
to  maintain  the  accuracy  for  the  numerical  integration  pro- 
cess, then  re-storing  every  other  point  of  the  1800  point 
arrays  as  900  point  arrays  to  accomodate  Subroutine  GRAPH. 
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4.     RESULTS  AND  DISCUSSION 

a.  Character  of  the  response 

An  understanding  of  the  results  is  dependent  upon  an 
appreciation  of  the  distinctive  features  of  the  loading  and 
their  effects  upon  the  response.   Loading  is  accomplished  by 
a  forced  displacement  of  one  end  of  the  column  in  the  form 
of  a  half-sine  pulse.   This  pulse  is  characterized  by  the 
dimensionless  parameters  c*-  and  G>  .   The  quantity  c^  is  the 
ratio  of  the  displacement  pulse  height  to  that  end  displace- 
ment which,  in  the  absence  of  lateral  deflection,  corresponds 
to  the  Suler  buckling  load.   The  parameter  G>  is  the  ratio  of 
the  first  mode  natural  period  of  lateral  vibration  of  the  col- 
umn model  to  the  pulse  duration.   Thus,  ex  measures  the  ampll. 
tude  of  the  loading  and  ^  measures  its  rapidity. 

Consider  the  very  slow  loading  ( Q> « I  )  of  a  slender 
column  to  a  maximum  end  displacement  beyond  the  Euler  criti- 
al  value  (o<>  I  )0   At  axial  loads  above  the  Euler  limit 
the  straight  form  of  the  column  i3  unstable*   To  Initiate 
lateral  deflection  of  the  column  model  the  end  pins  are 
eccentrically  located.   The  effect  of  this  eccentricity   is 
to  produce  bending  at  loads  below  the  stability  limit  of  the 
column.   The  maximum  force  (Pfn  )  developed  in  the  column 
becomes  less  than  the  Euler  load  (P^  ).   However,  for  very 
slow  loading  the  maximum  force  developed  in  the  column  model 
is  the  static  maximum  load.   Increasing  values  of  C<  above 
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that  corresponding  to  the  Euler  critical  value  do  not  change 
the  value  of  the  maximum  force  Pm . 

For  increasing  loading  rates  (increasing  values  of  0>  ) 
the  lateral  motion  of  the  column  is  retarded  by  the  inertia 
of  the  mass  and  any  meaningful  correspondence  to  static  be- 
havior is  lost.   However,  the  maximum  force  developed  in  the 
column  is  related  to  the  rate  of  loading.   Also,  under  pulse 
loading  the  inelastic  column  recovers  from  a  deflected  posi- 
tion with  an  amount  of  residual  deflection  and  retains  a 
capacity  to  sustain  further  load.   This  capacity  may  be 
measured  in  terms  of  the  amount  of  residual  deflection  y>- 
which  results  from  the  loading.  This  residual  deflection  is 
equivalent  to  additional  eccentricity  for  any  subsequent  load- 
ing.  The  maximum  force  Pm  and  the  residual  deflection  y»-  appear 
to  be  the  most  important  performance  criteria  fcr  the  dynamic 
behavior  of  the  inelastic  column  subjected  to  a  pulse  loading, 
Dimenslonless  forms  of  these  two  quantities  are  utilized  in 
the  graphical  form  of  the  results. 

b.  Range  of  variables 

Solutions  for  Pm/PE  (maximum  force)  and  for  yr/o   (residual 
deflection)  have  been  obtained  for  a  range  of  (>>  from  0.01  to 
1.0  and  for  a  range  of  ex*  dependent  on  the  column  slenderness 
ratio.   The  results  include  solutions  for  three  values  of  col- 
umn slenderness  ratio  (L/o  )  and  two  values  of  eccentricity 
ratio  (e/o  ).   For  slenderness  ratioa  of  70  and  50  the  mini- 
mum pulse  height  considered  was  c<=  lc   For  a  slenderness 
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ratio  of  100  the  minimum  pulse  height  considered  was  cx=   2. 
For  values  of  O^.  less  than  these,  inelastic  action  is  insig- 
nificant. The  maximum  values  of  ex.  considered  are  based  on 
the  maximum  value  of  tensile  stress  occurring  in  the  column 
model  springs.   For  the  investigation  the  stress- strain  re- 
lationship after  departure  from  the  compressive  stress-strain 
curve  is  considered  linear  into  the  tensile  stress  region. 
For  this  reason  the  maximum  tensile  stresses  in  the  column 
are  limited  to  values  near  the  nominal  yield  stress  of  the 
column  material.   For  slenderness  ratios  of  70  and  50  the  max- 
imum pulse  height  considered  was  c<  =  205.   For  a  slenderness 
ratio  of  100  the  maximum  pulse  height  considered  was  cx=*  4. 
For  values  of  ex  greater  than  these  the  maximum  value  of  ten- 
sile stress  occurring  in  the  column  model  is  considerably  in 
excess  of  the  nominal  yield  stress  for  some  rates  of  loading. 

c.  Maximum  force 

Figure  10  is  a  plot  of  families  of  curves  of  maximum  force 
Pm/Ps:  versus  G  for  a  constant  pulse  height  ex.  .   The  families 
include  slenderness  ratios  of  100,  70,  and  50  for  an  eccentri- 
city ratio  of  C.05.   Figure  11  is  a  plot  of  the  same  families 
of  curves  for  an  eccentricity  ratio  of  C.01.   The  lower  limit 
of  G>  =  0.01  is  for  slow  loading  rates  approaching  the  static 
case.   The  upper  limit  of  G)  =  1.0  is  approaching  the  region 
where  axial  Inertia  effects  become  important  in  the  elastic 
column  as  reported  by  Taylor  /5/«   As  the  Inelastic  column 
model  does  not  include  the  effects  of  axial  inertia  G>  =  1.0 
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is  the  upper  limit  considered  in  this  study. 

Figures  10  and  11  show  that  as  the  rate  of  loading  is 
increased  the  maximum  force  developed  in  the  column  increases. 
This  is  due  to  the  delay  in  the  development  of  lateral  deflec- 
tion because  of  lateral  inertia.   These  figures  also  indicate 
that  as  the  rate  of  loading  is  decreased  the  maximum  force  de- 
veloped in  the  column  approaches  a  limit  asymptotically.   In 
Fig.  10  this  limit  may  be  observed  for  each  slenderness  ratio 
as  all  of  the  curves  converge  to  a  single  value  at  6  =  0.01. 
This  limit  is  the  static  load  capacity  of  the  column  for  an 
eccentricity  ratio  of  0.05.   In  Fig.  11  this  limit  is  not 
attained.   The  spread  of  points  indicates  that  for  this  eccen- 
tricity (e/o  ■  0.01)  there  is  considerable  dynamic  effect  at 
0  =  0.01.   However,  the  curves  of  minimum  &<   for  each  slen- 
derness  ratio  at  Q>  =  0.01  appear  to  approach  a  limit  of  the 
static  load  capacity  of  the  column  fcr  an  eccentricity  ratio 
of  C.01.   A  comparison  of  these  figures  confirms  that  as  the 
eccentricity  increases  the  maximum  force  developed  in  the 
column  decreases. 

d„  Slow  loading 

It  is  possible  to  make  some  observations  concerning 
the  values  obtained  at  the  low  end  of  the  P>  range  by  con- 
sidering the  static  case  of  a  centrally  loaded  inelastic 
column.   For  this  ca3e  the  maximum  load  the  column  can  sup- 
port is  between  the  tangent  modulus  load  FT  and  the  reduced 
modulus  load  PR  .   In  terms  of  the  Euler  load  the  dlmensionless 
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tangent  modulus  load  is 

P      Ft 

JLl  =  J=Z  (23) 

v/here  ET  is   the   tangent   modulus  or  the   local    slope   of  the 
compressive    stress-strain  diagram.      The   tangent   modulus   de- 
fines  the  value  of   stress   at   which  bending  begins   for  the 
static    case   of   a   centrally   loaded    inelastic    column.      In 
terms   of   the  model    constants  developed   previously   this 
critical   stress   is 


CTcr   = 


_  TT4      ET 


8    (Ly^y- 

or   in  terras   of  diraensi onless    stress  this   becomes 

The   ratio   of  Young's  modulus  to  the  tangent   modulus   in  tenrs 
of   the    Ramberg-Osgood   parameters    is 

fT=1+M£o±L(J|_r  (25) 

where  n  =  10  (the  shape  parameter)  and  S-,  =  C,0C4  (the  dimen- 
sionless  compressive  yield  stress)  are  the  parameters  used 
in  the  investigation.   By  choosing  a  value  of  Scy-  and  sub- 
stituting this  for  S  in  equation  (25)  the  ratio  E/ET  is  de- 
termined.  By  substituting  the  values  of  this  ratio  and  Sa- 
in equation  (24)  the  slenderness  ratio  Is  determined.   In 
this  way  a  value  of  the  ratio  ET/E  ^s  determined  correspond- 
ing to  each  of  the  slenderness  ratios  100,  70,  and  50.   This 
value  of  £t/2  equals  PT  /Pc  from  equation  (23) 0 
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In  terms  of  the  Suler  load  the  dimensionless  reduced 
modulus  load  Is 

Pe    EL  (26) 

where  Er.  is  the  reduced  modulus.   The  reduced  modulus  de- 
fines a  maximum  value  of  stress  which  can  be  approached  but 
not  exceeded  for  the  static  case  of  a  centrally  loaded  inel- 
astic column.   In  terms  of  the  model  constants  developed  pre- 
viously this  critical  stress  is 

rr     -  "r1"4  E« 
cr-"  8  (L/^ 

or  in  terms  of  dimensionless  stress  this  becomes 

The  reduced  modulus  is  dependent  on  the  tangent  modulus  and 
the  geometry  of  the  column.   For  the  column  model  it  can 
be  shown  that  the  reduced  modulus  is 

r  -  a£Er_ 
Lr  "  FTTT 

or  in  dimensionless  form  this  becomes 

T-~TTZvT  <28> 

and  this  value  is  equal  to  Pr/Pj:  from  equation  (26).   By 
choosing  a  value  of  SCl-  and  substituting  this  for  3  in  equa- 
tion (25)  the  ratio  ET/E  Is  determined.   The  ratio  ET/E  also 
determines  S-r/E  from  equation  (28).   By  substituting  the  values 
of  Er/S  and  SCr  in  equation  (27)  the  slenderness  ratio  is  deter- 
mined.  In  this  way  a  value  of  the  ratio  ER/E  is  determined 
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corresponding  to  each  of  the  slenderness  ratios  100,  70,  and 
50.   This  value  of  ER/E  equals  PR/PE  from  equation  (26). 

The  values  of  the  tangent  modulus  load  and  the  reduced 
modulus  load  determined  from  these  two  methods  are: 

for  L/^  =  100,   PT/PE  =  1.000  and  PR/Pe  =  1.000; 

for  L/^  =•  70,   PT/Pt  =  0.969  and  pR /pt  =  C.977; 

for  L/<^  =  50,   PT/PE  =  0.649  and  PR/PF.  =  O.69I  . 
From  Figs.  10  and  11  it  Is  seen  that  in  most  cases  the  max- 
imum values  at  ^>  =  0.01  are  less  than  the  tangent  modulus 
load  although  the  spread  of  points  for  an  eccentricity  ratio 
of  0.01  indicates  that  dynamic  effects  are  more  pronounced 
the  smaller  the  eccentricity.   For  this  reason  some  of  tne 
maximum  values  are  greater  than  the  reduced  modulus  load. 

To  observe  further  the  effect  of  eccentricity  on  the 
load  capacity  of  the  column  at  very  slow  loading  rates, a 
study  was  made  to  determine  toe  behavior  of  the  column  under 
a  ramp  displacement  of  the  end  at  a  loading  rate  correspond- 
ing to  that  available  to  commercial  testing  machines.   In 
terms  of  the  parameters  previously  used  to  describe  the  pulse 
loading  of  the  column  the  ramp  displacement  of  the  head  is 

X  =  Jcx  Q>  T 

where  cx=  1.0  and  &  ■  0.01.   Computer  solution  results  in- 
cluding a  graphical  presentation  of  Strain  1  (  €,  )  and  Strain 
2  (  6Z  )  versus  force  (P/P£ )  were  obtained  for  the  three  column 
slenderness  ratios  of  100,  70,  and  50  for  eccentricity  ratios 
of  C.05,  0.01,  and  0.001.   The  graphical  results  shown  In  App.II 
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exhibit  a  general  resemblance  to  the  experimental  strain 
versus  force  carves  of  actual  column  test  data  displayed  by 
Shaniey  /l/.   From  the  results  given  In  Table  1  It  Is  seen 
that  for  eccentricity  ratios  of  C.05  and  0.01  the  maximum 
load  is  less  than  the  tangent  modulus  load  for  all  three 
slenderness  ratios  considered.   However,  the  maximum  load 
increases  with  a  decrease  in  the  eccentricity  ratio.   As 
the  eccentricity  ratio  is  decreased  to  0.001  the  maximum 
load  is  greater  than  the  tangent  modulus  load  but  less  than 
the  reduced  modulus  load  for  a  slenderness  ratio  of  50.   For 
a  slenderness  ratio  of  70  there  is  considerable  dynamic  effect 
and  the  maximum  force  exceeds  the  reduced  modulus  load  by  a 


Table  I  Computer  and  Theoretical  Static  Load  Comparison 
(For  Ramp  Loading  of  Column  Model) 


L/p^  =   100 

L/^  =   70 

L/^  =   50 

Static 
PT/PE 

1.000 

O.969 

O0649 
O.69I 

Static 
PR/PE 

1.000 

0.977 

Dynamic 
(Computer) 

P*/PE 

for 
e/o  =  0.05 

0.952 

0.839 

0.584 

for 
e/o  =   0.01 

0.983 

0.916 

0.633 

for 

e/P  =   0.001 

1.000 

O.98I 

C.665 
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small  amount  *   For  a  slenderness  ratio  of  100  there  Is  consider- 
able dynamic  effect  also  but  the  solution  oscillates  about  an 
average  maximum  value  of  P^/  Pe  =  1,0. 

The  tabulated  values  of  P^/Pe  for  an  eccentricity  ratio 
of  0.05  for  a  ramp  displacement  loading  are  in  agreement  with 
the  maximum  force  values  Pm/Pr  for  pulse  loading  at  ft  =  0.01 
in  Fig,  10.  The  strain  (  €.,  and  €.2,)  versus  force  (P/PE  )  curves 
for  the  ramp  displacement  loading  for  an  eccentricity  ratio 
of  0.05  shown  in  Appc  II  display  no  dynamic  behavior.   In  Fig. 
11  (e/o  =  0 .01),  where  the  dynamic  effects  are  more  pronounced 
at  G>  »  C.01  for  pulse  loadings  the  forces  on  the  <X=  1  curves 
for  slenderness  ratios  of  70  and  50  agree  with  these  tabulated 
results.   For  a  slenderness  ratio  of  100  the  smallest  value  of 
c*k  considered  for  pulse  loading  was  2  and  this  curve  appears 
to  approach  a  value  Pm/P£  ■  0.983  asymptotically „   The  strain 
versus  force  curves  for  the  ramp  displacement  loading  for  an 
eccentricity  ratio  of  0.01  shown  in  App„  II  display  an  insignif- 
icant amount  of  dynamic  behavior  which  does  not  affect  the 
values  of  maximum  force  Pm/PE  .   Therefore,  these  values  of 
P,- /P£  may  be  considered  the  static  load  capacity  of  the  inel- 
astic column  for  eccentricity  ratios  of  0.05  and  0.01. 

Appendix  III  contains  a  compilation  of  the  dynamic  re- 
sponse curves  stress  (Si  and  3.j_  )  versus  time  (T)  for  a  range 
of  £>  from  0.1  to  1,0  and  for  an  eccentricity  ratio  e/o  =  0.01. 
This  graphical  compilation  includes  the  solutions  for  the  min- 
imum and  maximum  values  of  pulse  height  considered  in  the 
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study  for  all  three  slenderness  ration  0  Included  with  the 
compilation  Is  a  discussion  of  the  more  prominent  features 
of  the  curves. 

e.  Residual  deflection 

Figures  12 9    13 ,  and  14  are  famlles  of  curves  showing 
residual  deflection  yr/p  versus  0>  for  a  constant  pulse 
height  Ov  and  for  e/o  -  0001  and  0oo5o   Figure  12  is  the 
family  for  L/C  =  100,  Figc  13  is  the  family  for  L/o  =  70, 
and  Flg0  14  is  the  family  for  L/o  =  50 0   The  points  shown 
on  these  curves  are  from  the  identical  computer  solution 
results  as  the  points  shown  in  Figsc  10  and  11 0 

From  Fig0  12  for  L/o  =  100  it  Is  shown  that  in  general 
the  values  of  residual  deflection  increase  with  the  rapidity 
of  loading  to  a  maximum  value  then  decrease  abruptly 0   This 
abrupt  decrease  is  because  of  the  delay  in  the  development  of 
lateral  deflection  due  to  lateral  inert ia0   From  Fig„  13 
for  L/Q^  =  70  and  Fig0  14  for  L/^  =  50  it  is  shown  that  in  gen= 
eral  the  values  of  residual  deflection  are  essentially  con- 
stant as  the  rapidity  of  loading  Increases  until  an  abrupt 
decrease  occurs  due  to  the  effect  of  lateral  inert ia0   The 
exception  for  ex  -   1  and  L/o  =  70  is  due  to  the  small  amount 
of  Inelastic  action  that  takes  place  for  this  pulse  height 
at  all  values  of  P)  »      The  amount  of  residual  deflection  may 
be  expressed  in  terms  of  the  amount  of  inelastic  action  that 
occurs,,   From  the  graphical  compilation  of  stress  versus  time 
curves  in  Appc  III  and  the  data  compiled  in  AppD  V  the  amount 


42 


0.5 


0.2 


0.1 


0.05 


7r/? 


0.02 


0.01 


0.005 


0.002 


0.001 


0.01    0.02 


0.05 


Fig.  12.  Values  of  residual  deflection  y*-/^  vs.  (^ 
(ex  constant,  e/c^ -0.01  and  0.05)  for  a 
column  slenderness  ratio  L/^  of  100 
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Values  of  residual  deflection  y>/p  vs»  ^ 
(cxconstant,    e/^»  0.01  and  0.05)    for  a 
column  slenderness  ratio  L/c^     of  70. 
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Fig.  14.   Values  of  residual  deflection  yv/£  vs. ^ 
(  ex.  constant,  e/<£=  0.01  and  0.05)  for  a 
column  slenderness  ratio  L/^  of  50. 


of  Inelastic  action  may  be  determined,,   For  a  slenderness 
ratio  of  100,  inelastic  action  in  the  column  Is  limited  to 
maximum  compressive  stress  levels  below  the  value  of  the 
yield  stress  for  all  solutions  considered  In  the  investiga- 
tion.  For  a  slenderness  ratio  of  50,  inelastic  action,  in 
general,  occurs  at  much  higher  maximum  compressive  stress 
levels-near  or  above  the  yield  stress0   For  a  slenderness 
ratio  of  70,  Inelastic  action  occurs  at  maximum  compressive 
stress  levels  intermediate  to  those  for  slenderness  ratios 
of  100  and  50, 

A.  comparison  of  Fig.  12,  13,  and  14  indicates  that  as 
the  slenderness  ratio  decreases  the  effect  of  eccentricity 
on  the  values  of  residual  deflection  is  less  pronounced. 
However,  decreasing  eccentricity  causes  a  shift  of  the  re- 
sidual deflection  curves  to  longer  pulse  durations  (smaller 
values  of  (3  )D   This  is  because  dynamic  behavior  is  more 
pronounced  the  smaller  the  eccentricity  and  lateral  inertia 
effects  become  more  important  at  slower  loading  rates.   It 
may  also  be  observed  from  these  residual  deflection  curves 
that  the  maximum  values  of  residual  deflection  do  not  necessa- 
rily occur  at  the  slowest  loading  rates 0   The  minimum  values 
of  residual  deflection  generally  occur  at  the  fastest  loading 
rates.   The  only  exceptions  are  for  the  columns  of  largest 
slenderness  ratio  (L/o  =  100)  and  largest  eccentricity 
(e/o  =  0.05).   However,  the  residual  deflections  appear  to 
follow  the  same  trend  but  are  limited  by  the  maximum  value  of 
the  parameter^  ■  1.0(  shortest  pulse  duration  considered  in  the 

investigation) * 
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5.    CONCLUSIONS 

From  the  results  discussed  In  the  preceding  section, 
the  following  conclusions  may  be  drawn  for  the  rates  of 
loading  considered  in  the  investigation: 

(a)  The  inelastic  column  model,  together  with  a  high 
speed  digital  computer  and  the  peripheral  plotting 
equipment,  provides  a  sound  and  reasonable  method 
for  the  study  of  the  dynamic  behavior  of  inelastic 
columns,, 

(b)  For  very  slow  loading  (long  pulse  duration)  the  model 
exhibits  behavior  in  agreement  with  Shanley ' s  postu- 
lateso 

(c)  As  eccentricity  Is  reduced,  the  lateral  response  is 
slowed 0  Accordingly,  columns  with  very  small  eccen- 
tricity show  significant  dynamic  effects  at  quite 
slow  loading  rates0 

(d)  Within  the  range  of  pulse  durations  studied,  the 
maximum  load  Increases  monotonically  with  decreas- 
ing pulse  duration, 

(e)  The  ratio  of  maximum  dynamic  load  to  static  load 
capacity  decreases  as  the  slenderness  ratio  is 
reduced 0 

(f)  Within  the  range  of  pulse  durations  studied,  the 
smallest  residual  deflections  generally  accompany 
the  shortest  pulse  durations0 
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APPENDIX  I 

Example  of  Graphical  Output  from  Prop-ram  COL  1 

Graphical  Notation 

The  computed  points  on  the  graphs  are  the  values  of  the 

dimension] ess  variables  developed  in  the  basic  equations. 

The  £raphical  notation  is  as  follows: 

Symbol!  on  Graph  Corresponding  Symbols  in 

the  basic  equations 

SLENDERNESS  L/^ 

ECC  e/^ 

ALPHA  ex. 

BETA  Q 


STRAIN  1 


v 


STRAIN  2  £t 

STRESS  1  S, 

STRESS  2  S- 

DEFLECTION  y/^ 

FORCE  P/PE 

TIME  T 

Each  of  the  curves  on  the  graph  Is  labeled  to  correspond 
to  the  terms  in  the  graph  title,,   The  x-axis  is  the  axis 
of  abscissas  and  the  y-axis  is  the  axis  of  ordinates.   The 
numbers  specifying  the  coordinates  are  Integer  numbers  locat- 
ed to  the  right  of  the  coordinate  on  the  x-axls  and  above  the 
coordinate  on  the  y-axis „   The  scale  quantity  is  the  value  of 
the  first  marked  coordinate  from  the  origin  in  either  the 
x-direction  or  the  y-direction  expressed  in  powers  of  ten 

1-1 


to  locate  the  decimal  point c   The  symbol  E  in  the  scale  value 
indicates  that  the  number  following  is  the  exponential  power 
of  ten  to  be  applied  to  the  basic  scale  value,  i.e.,  1.00E+C1 


is  1.00  x  10"1"' 


Discussion  of  the  curves 


On  pages  1-5  to  1-7  following  this  discussion  are  shown 
curves  of  £,  ,  £z  ,  S,  ,  Sz   ,  y/^  and  P/PE  versus  T0   These 
curves  are  the  solutions  for  the  Inelastic  column  model  dy- 
namic response  to  a  ram  pulse  loading  with  the  following  par- 
ameters: 

(a)  slenderness  ratio,  L/o  =  100 

(b)  eccentricity  ratio,  e/o  =  0o05 

(c)  ram  pulse  height, ex.  =  4.0 

(d)  ram  pulse  duration  In  terms  of  £>  =  0.2, 
Tf  -  2TT/^  -  10TT 

The  discussion  includes  some  of  the  more  prominent  features 
of  the  dynamic  behavior  of  the  column  model  under  ram  pulse 
loading. 

A  comparison  of  the  strain  versus  time  curves  with  the 
stress  versus  time  curves  indicates  the  amount  of  inelastic 
action  that  takes  place  in  Spring  1  and  Spring  2.   In  Spring  1 
the  compressive  stress  level  increases  above  that  in  Spring  2 
as  the  amount  of  bending  increases.   The  internal  bending 
moment  in  the  column  model  is  directly  proportional  to 
Sj  -  Sa  .   As  the  amount  of  bending  increases  Spring  2  de- 
parts from  the  compressive  stress-strain  curve  along  a  straight 
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line  with  unit  slope  and  continues  into  the  tensile  stress 
region.  The  values  of  stress  in  Spring  2  remain  on  this 
straight  line  for  the  remainder  of  the  solution,,   The  com- 
pressive stress  level  in  Spring  1  continues  to  Increase  to  a 
point  where  it  also  departs  from  the  compressive  stress-strair 
curve  along  a  straight  line  as  the  internal  bending  moment  in- 
creases and  the  force  decreases  to  zero»  This  force  is  dir- 
ectly proportional  to  S  +  S-c  When  the  force  becomes  zero 
the  end  of  the  column  has  departed  from  the  ram  which  exhibits 
a  half -sine  pulse  displacement  «   The  column  model  is  then  in 
a  state  of  free  vibration  until  the  end  again  makes  contact 
with  the  ramc   The  compressive  stress  level  in  Spring  1  again 
returns  to  values  on  the  compressive  stress-strain  curve 0 
This  alternate  departure  and  contact  between  the  column  and 
the  ram  continues  until  the  ram  displacement  pulse  nears  its 
termination,,  The  column  then  continues  in  a  state  of  free 
vibration  where  S«  -  -  S>0   The  amount  of  residual  deflec- 
tion may  then  be  determined  from  the  mean  value  of  the  deflec~ 
tion  or  from  the  deflection  at  the  time  when  S>  *  S2  ■  0„   At 
this  point  the  internal  forces  and  the  internal  bending  moment 
in  the  column  model  are  zeroc 

The  effect  of  eccentricity  on  the  dynamic  behavior  of 
the  column  model  with  all  other  parameters  remaining  the 
same  may  be  observed  by  comparing  the  stress  versus  time 
curves  in  this  appendix  with  the  stress  versus  time  curves 
on  page  111-13  where  the  eccentricity  ratio  is  0o01o   It  may 
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be  noted  Miat  the  essential  differences  in  these  curves  are 
that  the  maximum  compressive  stress  levels  attained  in  Spring  1 
and  Spring  2  are  greater  for  the  smaller  eccentricity „   Also, 
contact  between  the  ram  and  the  column  terminates  sooner  for 
the  smaller  eccentricity,. 
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APPENDIX  II 

Theoretical  and  Computer  Solution  Graphical  Comparisons  of 
Maximum  Loads  for  the  Static  Case 

The  notation  for  the  graphs  displayed  in  this  appendix 
Is  explained  in  App0  I0 

The  graphs  on  pages  II<=3  to  11=11  folowlng  this  dis- 
cussion are  the  curves  I ,  and   ,  versus  P/F  for  the  three 
slenderness  ratios  considered  in  the  study 0   As  previously 
mentioned ,  the  loading  for  the  study  in  this  appendix  is  a 
ramp  displacement  of  the  end  of  the  column  model  at  a  dis= 
placement  rate  comparable  to  that  in  commercial  testing  ma= 
chines o   The  study  includes  eccentricity  ratios  of  0o05, 
OoOl,  and  0oC01  for  each  of  the  slenderness  ratios 0  The 
tabulated  results  of  this  study  are  shown  in  Table  I  in 
Section  40 

In  addition  to  the  previous  comments  made  concerning 
the  result b9    it  may  also  be  seen  from  the  graphs  that  for 
all  slenderness  and  eccentricity  ratios  considered  an  amount 
of  bending  takes  place  before  strain  reversal  occurs  in 
Spring  2  and  before  the  maximum  load  is  reached c   This  was 
essentially  Shanley°s  point  in  proving  the  fallacy  of  the 
reduced- modulus  theory  for  column  actionc   As  stated  by 
Shanley  /!/ s 


O     O  O  O  0 


othe  reduced=modulus  theory  is  not  correct 
for  predicting  the  maximum  load  up  to  which  a  per= 
feet  column  will  remain  straight 0   This  is  because 
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it  is  possible  for  the  column  to  bend  simultaneously 
Kith  increasing  axial  loadG   Under  such  conditions 
it  Is  possible  to  have  bending  without  introducing 
any  strain  reversal ,  upon  which  the  reduced- modulus 
theory  depends 0 
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APPENDIX  III 
Column  Model  Dynamic  Response  Curves 

The  notation  for  the  graphs  displayed  in  this  appendix 
is  explained  in  App.  I. 

The  graphical  compilation  on  pages  III-3  to  111-26 
following  this  discussion  gives  the  curves  S|  and  Sz  versus 
T  for  ©  =  0.1,  0.2,  0.5  and  1.0  and  for  an  eccentricity  ratio 
e/o  =  0.01.   Included  in  the  compilation  are  the  solutions 
for  the  minimum  and  maximum  values  of  ram  pulse  height  ex.  con- 
sidered in  the  investigation  for  all  three  slenderness  ratios. 
The  graphical  solutions  are  grouped  by  slenderness  ratio  with 
the  solution  for  the  maximum  value  of  c<  followed  by  the  solu- 
tion for  the  minimum  value  of  oc  for  a  single  &  value.   The 
order  of  the  &  values  is  0.1,  0.2,  0.5  and  1.0,  that  is,  from 
the  longest  pulse  duration  of  loading  to  the  shortest  pulse 
duration  of  loading  considered  in  the  compilation. 

A  discussion  of  some  of  the  prominent  features  of  the 
dynamic  behavior  of  the  column  model  under  ram  pulse  loading 
is  Included  in  App.  I  where  the  curves  of  6,  ,  6L,  S,  ,  SZi  y/o 
and  P/PE  versus  T  are  displayed  for  one  computer  solution. 
The  following  is  a  discussion  of  some  of  the  more  prominent 
features  of  the  dynamic  behavior  of  the  column  model  in  terms 
of  the  stress  curve  compilation  in  this  appendix. 

In  general  the  amount  of  inelastic  action  in  the  column 
model  is  inversely  related  to  the  slenderness  ratio.   For  a 
slenderness  ratio  of  100  Inelastic  action  is  limited  to  com- 
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pressive  stress  lsvels  below  the  value  of  the  yield  stress. 
For  a  slendemess  ratio  of  50  the  maximum  conpressive  stress 
levels  are  near  or  above  the  yield  stress.   For  a  slenderness 
ratio  of  70  inelastic  action  occurs  at  maximum  compressive 
stress  levels  intermediate  to  those  for  slenderness  ratios 
of  50  and  100. 

It  may  also  be  noted  that,  in  general,  contact  is  lost 
between  the  ram  and  the  end  of  the  column  (when  S,  =  -  SL) 
before  the  loading  pulse  nears  its  termination  for  ft  values 
of  0.2  and  0.5  and  slenderness  ratios  of  70  and  50.   For  a 
slenderness  ratio  of  100  this  occurs  at  values  of  C>  =0.1, 
0.2  and  0.5. 

It  may  be  seen  that  the  value  of  maximum  tensile  stress 
in  Spring  2  exceeded  the  value  of  the  nominal  yield  stress 
(S^j  -  C.CC4)  for  the  C>  values  of  C.2  and  0.5  and  ex.  -  2.5 
for  a  slenderness  ratio  of  50. 
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APPENDIX  _IV 

Computer  Program s  and  Not at  Ion 

a.    General 

The  FORTRAN  language  computer  programs,  together  with 
the  notation,  and  the  program  operating  Instructions  for 
both  Program  COL  1  and  Program  COL  2  are  Included  In  the 
following  pages.   The  discussion  and  a  listing  of  Program 
COL  1  are  included  on  pages  IV  -  2  to  IV  -10.   The  discuss- 
ion and  a  listing  of  Program  COL  2  are  included  on  pages 
IV  -11  to  IV  -20.   The  FORTRAN  language  format  for  each  of 
the  programs  is  in  standard  FORT, RAN  language  to  be  used  with 
a  standard  FORTRAN  compiling  routine,   The  program  language 
is  compatible  with  any  of  the  large  machine  FORTRAN  compiling 
systems  including  the  I3M  7090  and  the  CDC  16C4  systems „   Each 
line  in  the  listing  is  the  information  contained  on  one  80  col- 
umn IBM  card  where  columns  1-5  contain  any  statement  number, 
column  6  is  used  to  indicate  a  continuation  card,  columns 
7-7?  contain  the  FORTRAN  statement,  and  columns  73-80  are 
used  for  card  identification  only.  Although  the  notation  is 
basically  the  same,  the  subscript  notation  is  different  for 
each  program.   Program  COL  1  stores  a  number  of  varlabls 
quantities  in  array  form  for  print-out  and  to  accommodate 
Subroutine  GRAPH,  whereas,  Program  CCL  2  stores  only  one 
value  of  tne  variable  quantities  for  print-out. 


IV-1 


b„  Notation  for  Program  COL  1 

In  the  notation  for  Program  COL  1  the  subscript  J  re- 
fers to  quantities  within  the  Runge-Kutta  prediction  loop 
which  are  calculated  for  intermediate  times  within  the  time 
increment  AT.   The  subscript  I  refers  to  the  quantities 
corresponding  to  each  time  T(I).   The  notation  for  Program 
COL  1  is  as  follows: 

Computer  program  notation  for  Program  COL JL 


Computer 
Program 
_Term_ 

A 

ALPHA 

BETA 

DDV(J) 

DEFL(I) 

DELT 

DISPL(I) 

ECC 

EPSIl(J) 

EPSI2(J) 

FORCE(I) 


Corresponding 
terms  in  table 

_qH-3X^oli§ 

n 
ex. 

V 

V 

AT 
x/xE 
e/<^» 

P/PE 


Computer 
Program 
term  _ 

ST^xY 

SIG1(J) 

SIG2(J) 

SLNESS 

STRNl(I) 

STRN2(I) 

STRSl(I) 

3TRS2(I) 

T(I) 

VEL 
X(J) 


Corresponding 
terms  in  table 
of  symbols _ 


3; 

L/<^ 

€, 

s, 

Si 

r 
v 

x/xE 


c.    Operating   instructions   for  Program  COL  1 

The  output    from  Program  COL  1    includes   the  print-out 
and   the   array   quantities   necessary   for  Subroutine   GRAPH. 
The  listing  for  Program  CCL  1  does  not    include   Subroutine 
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GRAPH.      However,    a  listing  of   Subroutine   GRAPH  and   instruc- 
tions  for  its   use  may   be   obtained   from  the   Control   Data   Cor- 
poration.     To   use  the  program  to   obtain  the  print-out   quanti- 
ties  without   the   graphical   output    remove  the   cards   CALL  GRAPH. 
To  operate  Program  COL  1  without   Subroutine   GRAPH  use  the   cards 
including  Subroutine   STRESS   in  the   order   shown    in  the   listing 
followed   by   one   blank   card,    then  the    input   data  card. 

The   read-write   statements   in   the  listing  are  written 
for  magnetic   tape   input   and   output   as   follows: 

(1)  Input    (READ   INPUT   TAPE   2) 

The  input  to  the  program  Is  contained  on  one  card  (7E1C.3) 

and  includes  the  following  seven  quantities  in  order: 

Input  quantity  Format 

SLNESS  (L/<^  )  E10.3 

ALPHA  (  ex  )  E10.3 

ECC  (e/^)  E10.3 

BETA  (  qb  )  E10o3 

A   (n)  E10.3 

SIGY  (3b)  ElCo3 

DELT  (AT)  E10„3 

If  Subroutine  GRAPH  is  used  axis  scaling  and  multiple 

curve  data,  graph  title  data  and  curve  label  data  in 

addition  to  the  above  are  also  needed  in  the  input. 

(2)  Output    (WRITE  OUTPUT   TAPE  4) 

The  output  includes  the  input  quantities  and  the  vari- 
able quantities  for  every  fifteenth  point  as  follow: 
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Output   Quantity  Format 

SLNSSS    (L/^  )  F5»0 

ALPHA    (  ox.  )  F5o2 

ECC    (e/^>  )  F5o3 

BETA    (  Q>   )  24.2 

A      (n)  25*1 

3IGY    (Sb)  F603 

DSLT    (AT)  F5o3 

T  P6.2 

DISPL    (x/xt)  214 04 

STRN1    (£,)  214  „  4 

STRN2    (£..)  214.4 

ST  RSI    (S,  )  214  04 

3TRS2    (Sj  214 04 

DEFL    (y/^  )  214  c4 
FORCE    (P/Pr)  214 04 


i. 
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PROGRAM  COL  1 
C  PROGRAM  TC  DETERMINE  INELASTIC  COLUMN  RESPONSE  TO  HALF-SINE 
C  RAM  DISPLACEMENT  PULSE  APPLIEC  AT  END  WITH  SMALL  ECCENTRICITY 
C  TO  INITIALLY  STRAIGHT  COLUMN.  WHEN  RAM  CAUSES  TENSION  COLUMN 
C  EXHIBITS  FREE  VIBRATION. 

ODI MENS  ION  T(900), CI SPL ( 9CC ) , CEFL ( 900 ) ,STRN 1 C9C0) ,STRN2(900), 
1STRS1  (90C),STRS2(  90C)  ,  X(  4  )  ,  QCV<  U  )  f  EPS  1 1 C  ** )  ,  EPS  12  (U  *  •  S  IG  1  ( «i ) . 
2SIG2(U),AK(U,2),C ( U  ), FORCE (900) 

COMMON  A,SIGY 

READ  INPUT  TAPE  2 , 5 ,SLNESS , ALPHA, ECCt BETA, A, S IGY.DELT 
5  F0RMAT(7E1C3) 

T( 1)=C.C 

DEFLI 1)=0.0 

VEL=O.C 

STRN1 (0)=0.0 

STRN21C)=0.0 

STRS1  (0=0.0 

STRS2(C)=0.C 

C<1 )=C.O 

C(2)*C.S 

C(3>  =  C5 

C(U)=1.C 

KTENS=1 

K0NTACT=1 

M=0 
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1  =  1 

UO  DO  70  J«1,U 

TC=T(  I)4C( J)»DELT 

VC=DEFL(  I  )+C( J)»AK< J- 1 , 1 ) 

IF(0.5*BETA«TC-3.  1U 1 593) 50, 55, 55 
50  XU)  =  />LPHA#SINF(0.5«BETA«TC) 

GO  TO  60 
55  X(J)=C.O 

600EPSI1 (J)=12.1761/ (  SLNESS»«2  )  •< X ( J)-0. 5* ( ECC««2 )• VC««2* 
1E-COVCM  1.0-ECO) 

IF( J-1)62,62,63 

62  CALL  STRESSlEPSil U),STRN1(I-1),STRS1 (I-1),SIG1( J)  ) 
STRN1  (I  )=EPSI 1 (1) 

STRS1  ( I  )=SIG1  (1) 
GO  TO  6k 

63  CALL  STRESS(EPSI1 ( J) ,STRN1< I ),STRS1 <I),SIG1(J)) 
6U0EPSI2(J)=12.1761/ (  SLNESS»*2  ) • ( X( J )-0.5» ( ECC««2 )♦ VC«»2- 

1ECC»VC*(  1.0+ECC)) 
IF( J-1)65,65,66 

65  CALL  STRESS(EPSI2( J),STRN2( 1-1) ,STRS2( I- 1 ),SIG2( J)  ) 
STRN2( I  )  =  EPSI2(1) 

STRS2( I  )  =  SIG2(1 ) 
GO  TO  661 

66  CALL  STRESS(EPSI2(J) , STRN2( I )  ,  STRS2< I ),SIG2(J)) 
661  IF(KTENS-1  )67,665,665 
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665  IF(K0NTACT-l)67,69Ut69U 

67  IF(M-1)675,68,685 
675  X1=X(J) 

P1*SIG1CJ)+SIG2(J) 
X(J)=CISPL( 1-1) 


GO  TO  60 
68  X2-XU) 

P2=SIG1( J)*SIG2(J) 

X { J  )  =X2- i P2 / { P2-P 1 ) ) • ( X2-X 1 ) 

GO  TO  6C 

685  P=SIG1(  J)+SIG2(J> 
IF(P)686,688,688 

686  IF( 10.0»* ( -8. O)-AeSF(PJ) 6  87,693,693 

687  X1=X(J) 
P1=P 

X( J)=X2-(P2/(P2-Pll)»(X2-Xl) 
GO  TO  60 

688  IF(10.0«»(-8.0)-ABSF(P))689,693,693 

689  X2=X(J) 
P2=P 

X(  J)  =  X2-(P2/( P2-P1) )«(X2-X1) 
GO  TO  6C 
693  *»0 
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KTENS=1 
GO  TO  698 
69U  IFCSIGK  J)*SIG2(J )) 695 ,698, 698 

695  IF(C.5«eETA«TC-3. 1U 1593)6955, 696,696 
6955  KTENS=0 

GO  TO  675 

696  KONTACT=C 
GO  TO  675 

6980DDV(J)=C.0U10639»<    SLNESS«»2    ) •< ( SIG1 ( J)+SIG2 ( J ) )»(VC*1.0)- 
1 (S IG1  (J)-SIG2U))/ECC) 
AK(  J,  l)=DELT*(VEL-»-C<J)«DELT»COVU)  ) 
70    AK( J,2)=CELT»D0V(  J) 
DISPLt I )=X(1) 

FORCEl  I  )=.0410639M    SLNESS»«2     ) •( STRS 1 ( I )*STR$2 (  I )  ) 
VEL  =  VEL*(AK(1,2)*2.0»AK(2,2)*2.0«AK(3,2)*AM4,2)  1/6.0 
DEFL(  14 1  )=DEFL( I)  ♦< AK ( 1 , 1 )*2.0«AK< 2, 1 )*2.0«AK (3, 1 ) *AK { U, 1  )  )/6.0 
IF( 1-900)75,80,80 
75    T(I  +  1  )=T(I)*DELT 

1=1*1 

GO  TO  UC 
80  WRITE  OUTPUT  TAPE  U,82 
82  F0RMAT11H1) 

WRITE  OUTPUT  TAPE  U,85,SLNESS, ALPHA, ECC, BETA, A, S IG Y.CELT 
850F0RMAT(  1UH    OUTPUT  DATA  //3X 12HSLENDERNESS«F5.0,2X6HALPHA=F5 .2, 

12XUHECC=F5.3,2X5H8ETA=FU.2,2X2HA*F5.1,2X5HSiGY«F6. 3,2X5H0ELT*F5. 3 

2    //)  rt 
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WRITE  OUTPUT  TAPE  U,90 
9CCF0RMAT(14XUHTIME,«*X12HDISPLACEMENT,UX10H0EFLECTIGN,UX8HSTRAIN  1  , 
16X8HSTRESS  1.6X8HSTRAIN  2.6X8HSTRESS  2.6X5HFORCE// ) 

DO  92  1=1,900 
92  DEFL( I)=6CC»DEFL<  I) 

DO  95  I=1,9C0,15 
950WRITE  OUTPUT  TAPE  U , 98, T ( I ) , D ISPLC I) , DEFL ( I ) , STRM ( I  )  , STRS  1  <  I )  , 

1STRN2U  ),STRS2(I)  ,FCRCE(  I) 
98  F0RPAT(3XF6.2,7E1U.U) 

CALL  GRAPH(900,T,  STRM, 8) 

CALL  GRAPH(900,T, STRN2,8) 

CALL  GRAPH(9C0,T, STRS1.8) 

CALL  GRAPH(900,T, STRS2.8) 

CALL  GRAPH(900,T,DEFL,8) 

CALL  GRAPH(900,T,  FORCE, 8) 

STOP 

END 

* 

SUBROLTINE  STRESS ( EPS  I ♦ EPS  I  A, SIGA.S IG ) 
COMMON  A.SIGY 
I FCEPS I )20, 20,21 

20  SIG=SIGA*(EPSI-EPSIA) 
GO    TO    33 

21  IF(EPSI-EPSIA)22,22,23 

22  SIG*SIGA-MEPSI-EPSIA> 
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GO    TO    33 

i 

23    SIGB=SIGA*(EPSI-EPSIA) 

B  =  3.0/(7.0«(SIGY»«(A-1.0n) 

SIG1=EPSI 
2U    SIG2=SIG1-<SIGUB»(SIG1«*A)-EPSI)/(A»B»(SIG1»»(A-1.0> )*1.0) 

IF(ABSF(SIG2M10. 0»»(-6.0>))-ABSF(SIG2-SIGl )) 29, 30 ,30 

29  SIG1=SIG2 

GO    TO    21*  .. 

30  IF(SIG2-SIG8)31.31,32 

31  SIG=SIG2 
GO    TO    33 

32  SIG=SIGB 

33  RETURN 
END 
ENO 
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d.  Notation  for  Program  COL  2 

In  the  notation  for  Program  COL  2  the  subscript  J  re- 
fers to  quantities  within  the  Runge-Kutta  prediction  loop 
which  are  calculated  for  intermediate  times  within  the  time 
increment  AT.   The  variable  quantities  ending  in  B  corre- 
spond to  a  time  T  -  AT  and  the  variable  quantities  ending 
in  A  correspond  to  a  time  T.   The  notation  for  Program  COL  2 
is  as  follows: 

Computer  program  notation  for  Program  COL  2 


Computer 
Program 
Term 

Corresponding 
terms  in  table 
of  symbols 

Computer 
Program 
Term 

Correspond lng 
terms  in  table 
of  symbols 

A 

n 

SIG1(J) 

3, 

ALPHA 

ex 

SIG2(J) 

S£ 

BETA 

Q> 

SLNESS 

L/^ 

DDV(J) 

«  » 

V 

STRN1A  or 

B 

€. 

DEFLA  or  B 

V 

STRN2A  or 

3 

€* 

DEFLP 

y*-/^ 

ST RSI A  or 

B 

s, 

DELT 

AT 

STRS2A  or 

3 

s* 

DISPLA  or  B 

sec 

x/x£ 
e/cD 

STRS1M 
STRS2M 

Max 
or  i 

Max 
or  < 

S,  (tensile 
compressive) 

3Z    (tensile 
compressive 

EPSIl(J) 

€, 

T 

T 

EPSI2(J) 

€* 

VELA  or  B 

V 

PMAX 

P^/Pe 

x(J) 

x/x£ 

SIGY 

SM 
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e.   Operating  Instructions  for    Program  COL  2 

To  operate  Program  COL  2  use  the  cards   including  Sub- 
routine  STRESS  in  the  order  shown   in  the  listing  followed   by 
one  blank  card,   then  the   input  data  card. 

The   read-write   statements  in  the  listing  are  written  for 
magnetic   tape  input   and  output   as  follows: 

(1)  Input  (READ  INPUT  TAPS  2) 

The  input  to  the  program  is  contained  on  one  card  (7E10.3) 

and  includes  the  following  seven  quantities  in  order: 

I nput  quant  1 1 y  Format 

SLNESS(L/^  )  S10.3 

ALPHA  (ex.  )  E10.3 

ECC  (e/cO  E10.3 

BETA  (  Q   )  E10.3 

A  (n)  E10.3 

SIGY  (Sa)  E10.3 

DELT  (AT)  E10.3 

(2)  Output  (WRITE  OUTPUT  TAPE  4) 

The  output  includes  the  input  quatities  and  one  value 
of  each  of  the  variable  quantities  Including  the  time 
T  of  occurrence  of  the  value  as  follows: 
Output  quantity  Format 

SLNE3S  (L/o)  F5.0 

ALPHA  (ex.  )  E14.4 

ECC  (e/^  )  F5.2 

BETA  (  Q  )  F4.2 

A   (n)  F5.1 
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Output   quantity  Format 

SIGY   (Sy)  F6.3 

DELT    (AT)  F5.3 


PMAX    (Prn/Pj.)  E14.4 

Tl  F6.2 


STRS1M   (Slmax)  E14.4 

T2  F6.2 


STRS2M   (SZmax)  E14.4 

T3  F6.2 


DEFLP    (Vr/p)  E14.4 

T4  n  F6.2 
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PROGRAM  COL  2 
C  PROGRAM  TO  DETERMINE  RESIDUAL  DEFLECTION  AND  MAXIMUM  FORCE 
C  IN  INELASTIC  COLUMN  RESULTING  FRCM  HALF-SINE  RAM  DISPLACEMENT 
C  PULSE  APPLIED  AT  END  WITH  SMALL  ECCENTRICITY  TC  INITIALLY 
C  STRAIGHT  COLUMN.  WHEN  RAM  CAUSES  TENSION  COLUMN  EXHIBITS  FREE 
C  VIBRATION. 

OD I  MENS  ION  X(W),DDV(lt)vEPSI1(4)vEPSI2(4)vSIGM«)«SIG2(U).AK(U,2)  , 
1C(U) 

COMMON  A.SIGY 

READ  INPUT  TAPE  2 ,5,SLNESS, ALPHA, ECC. BET  A, A,S ICY,DELT 
5  FORMAT(7E10.3) 
10  T=0.0 

DISPLA=0.0 

VELA=0.0 

OEFLA=C.C 

STRNie=0.0 

STRN2e=0.0 

STRS1E=C.C 

STRS2e=0.0 

PMAX=0.0 

STRS1M=C.C 

STRS2M=C.O 

STRS2E-0.0 

WORK=C.O 

C(1)=C.C 
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C(2)=C5 

C13)=0.5 

C(U)»1.0 

KENERG»1 

KTENS^l 

K0NTACT=1 

M=0 
UO    DO    70    J=l,4 

TC=T*C(J)*DELT 

VC-DEFLA*CU)«AK< J-1,  1) 

IF(0.5»BETA»TC-3,  114  1593)50,55,55 
50    X( J)=ALPHA*SINF(0.5«BETA»TC) 

GO    TO    60 
55    X(J)=C.C 

600EPSI1  (J)  =  12.1761/ (    SLNESS»»2    )  •(  X(  J  )-0.  5«(  ECC«*2  )•  VO«2* 
1ECC«VC*(1.0-ECC)) 

IF( J-l)62,62,63 

62  CALL    STRESS(EPSI1  ( J ) , STRN1B, STRS 1B,SI G1 (J) ) 
STRN1A*EPSI1(  1) 

STRS1A=SIG1(1 ) 
GO    TO    64 

63  CALL    STRESS(EPSI1 (  J  )  ,  STRN1  A,  STRS  1  A,SI  Gl  U)  ) 
640EPSI2(J)=12.1761/<    SLNESS»»2    )•  ( XU  )-0.5»(ECC»«2)»  VC«»2- 

.   1ECC*VC«<1«0+ECC)) 
IF (J- 1)65,65,66 
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65  CALL    STRESS(EPSI2(J),STRN2B,STRS2B,SIG2< J)) 
STRN2/>  =  EPSI2(1) 

STRS2A=SIG2<1) 
GO    TO    661 

66  CALL    STRESS(EPSI2(J),STRN2A, STRS2A.SIG2C J)) 
661     IF(KTENS-1 )67,665,665 

665     IF(KOKTACT-l)67,69U,69U 

67  IF(M-1)675, 68,685 
675    X1=X(J) 

P1=SIG1(J)  +  SIG2U) 

X(J)=CISPL8 

M  =  M  +  1 

GO    TO    6C 

68  X2=X(J) 

P2  =  SIG1 ( J)*SIG2(J  ) 

X( J)=X2-(P2/( P2-P1) >«<X2-X1) 

M=M*1 

GO    TO    60 

685  P=SIG1 ( J)+SIG2«J> 
IF(P)686,688,668 

686  IF( 10.C** (-8.0 )-ABSF«P 1)687,693,693 

687  X1=X(J) 
P1=P 

X(  J)=>2-(P2/(P2-P1))«(X2-X1 J 
GO    TO    60 
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688  IF  (10.0«M-8.0)-ABSF(P)  1689,693,693 

689  X2=X(J) 
P2  =  P 

X( J)=X2-(P2/(P2-P1))«(X2-X1) 

GO  TO  6C 
693  N=0 

KTENS=1 

GO  TO  698 
69U  IF(SIG1(J)+SIG2(J ) )695,698, 698 

695  IF(0.5»BETA»TC-3.  1U  1593)  6955,  696,696 
6955  KTENS=0 

GO  TO  675 

696  KONTACT=C 
GO  TO  675 

6980DDV< J)=0.0U10639»  (  SLNESS»«2  )«USIG1  ( J  )+S<It»2(  J)  )»CVO1.0- 
1 (SIG1 (J)-SIG2(J))/ECC) 
AK( J, 1)=CELT«(VELA+C( J) »DELT«DDV( J ) ) 

70  AK( J,2)=DELT»DDV( J) 
IF(0.5«BETA«T-3.1 Ul 593 )7C2,70U, 70U 

7020V*ORK  =  WORK  +  0.5M<STRS1A*STRS2A)*(STRS1B*STRS2B))» 

1 (ALPHA*(SINF(0.5»BETA»T)-SINF{0.5«BETA«(T-DELT>)  )) 
704    IF(PM/>X-C.0410639»<    SLNESS»»2    )  •  (STRS1  A+STRS2A  )  )71  ,72  ,72 

71  PMAX=C.CU10639M    SLNESS»»2    ) * ( STRS1 A+STRS2A) 
T1  =  T 

72  IF(STPS1V-STRS1A)73,7U,7U 
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73    STRS1MSTRS1A 

STRS1E*STRS1M 

T2  =  T 
7U    IF(A8SF<STRS2M)-ABSF{STRS2A>)75,751,751 

75  STRS2N=STRS2A 
T3  =  T 

751  IF(STRS2E-STRS2A> 752,755,755 

752  STRS2E=STRS2A 

755  IF(KENERG-1)756, 76,76 

756  IF(VELB)758,759,7'57 

757  IF(VELA)759,759,8.1 

758  IF(VELA)81,759,759 

759  B  =  3.0/(7.0*(SIGY»MA-1.0))) 

ENERG1=  (A»B/(A+1  .0)  )  •STRS  1  £••  ( A*  1  .  0)  ♦C.5»SrTRS  1B«»2 
ENERG2=  (A*8/(A*1 .0) ) «STRS2E»* ( A*l .0) +0. 5*STRS2E««2 
ENERGT=0. 0821278* (  SLNESS«»2  >  •<  ENERGUENERG2 ) 
GO  TO  eu 

76  IF(K0NTACT-1)77,81,81 

77  BEN0M1=STRS1B-STRS2B 
BEN0M2=STRS1 A-STRS2A 
IF(STRS1B)78,78,79 

78  IF(STRS1A)81,80,80 

79  IF(STRS1A)80,80,81 

80  DEFLP=(DEFLA    -BEN CN2/ ( BENDM2-BENDM 1 ) • ( DEFLA-OEFLB) MECC 
TU  =  T 
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KENERG=0 
81  T*T+DELT 

DISPLe-OISPLA 

DISPLA=X<  1) 

VELB=VELA 

VELA=\/ELe  +  (flK(i  ,2  )+2.0«AM2,2)  +  2.0*AK(3,2)*AKU,2)  )/6.0 

OEFLB=DEFLA 

OEFLA   ■   OEFLB   ♦( AKC 1 , 1  )  ♦2.0«AKC 2, 1  )*2. 0»AK (3,  1) *AK ( U  ,  1  )  J/6  .  C 

STRN1B=STRN1A 

STRN2E=STRN2A 

STRS1E=STRS1A 

STRS2B=STRS2A 

GO    TO    UO 
8k    WRITE    OUTPUT    TAPE    U,85,SLNESS,    ECC    .BETA , A , $IGY,OELT 
850FORMAT(///1UH   OUTPUT    OATA       //12HSLENDERNESS«F5.0,2X6H      ECOF5.2, 
12X5HBETA=FU.2,2X2HA=F5. l,2*5HSIGY*F6.3,2X5HOELT»F5.3//) 

WRITE    OUTPUT    TAPE    Uf86 
860FORNAT(6X5HALPHA, 7X1 1HOEFORMAT ION, 3X9HMAX    FORCE, 5X12HMAX    STRESS    1, 
12X12HNAX    STRESS   2//) 

87  WRITE    OUTPUT    TAPE    4,88, ALPHA, OEFLP,PMAX, STRSlMrSTRS2M 

88  F0RMAT(5E1M.U) 

WRITE  OUTPUT  TAPE  4, 885, TU, Tl , T2, T3  . 
885  FORMAT ( 20X 1H( F6.2 , 1H) ,6X1 H(F6 . 2, 1H) ,6XlH (F6.2, 1H ), 6X1 H(F6.2,  1H)  ) 

WRITE  OUTPUT  TAPE  4 ,887, WCRK, ENERGT 
887  FORMAT(25X5HWORK*E12.4,lCX7HENERGY«E12.4/) 
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STOP 

ENO 

SUBROUTINE  STRESS < EPS  1 1 EPSIA ,SIGA,S 16 ) 

COMMON  A.SIGY 

IF(EPSI>20,20,21 

20  SIG=SIGA+<EPSI-EPSIA) 
GO    TO    33 

21  IF(EPSI-EPSIA)22,22,23 

22  SIG=SIGA*(EPSI-EPSIA) 
GO    TO    33 

23  SIGB^SIGA-MEPSI-EPSIA) 
B=3.0/(7.0«(SIGY»MA-1.0))) 
SIG1*EPSI 

2U    SIG2  =  SIGl-(SIGl*B«(SIGl««A)-EPSi)/CA«B»(SIG1»«MA-l.C)MK0) 
IF(ABSFtSIG2«(  10.  C««(-8.Cm-ABSF(SIG2-SIGU  J  29,30,  30 

29  SIG1=SIG2 
GO    TO    2U 

30  IF(SIG2-SIGB)31, 31,32 

31  SIG*SIG2 
GO  TO  33 

32  SIG=SIGB 

33  RETURN 

END 

... 

END 
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APPENDIX  V 
Maximum  Force.  Residual  Deflection  and  Maximum  Stress  Data 

Tabulated  on  the  following  pages  are  the  maximum  force, 
residual  deflection,  and  maximum  stress  data  for  the  three 
slenderness  ratios  and  two  eccentricity  ratios  for  all  values 
of  0«*  and  G  which  were  Investigated  for  pulse  loading*  The 
values  of  maximum  lateral  deflection  that  were  available  are 
also  Included  in  the  tabulation.  These  values  were  only  avail- 
able from  the  solutions  to  computer  Program  COL  1  used  to  ob- 
tain graphical  results.  Thus,  only  a  limited  number  of  these 
values  are  included  in  the  tabulation. 

The  dimensionless  time  Increments  AT  used  for  the  numer- 
ical integration  were  the  same  for  a  given  value  of  G)  in 
all  the  solutions  for  both  computer  programs.  These  values 
are  as  follows; 

ft  AT 

0.01  0.10 

0.05  0.10 

0.10  0.08 

0.20  0.04 

0.50  0.02 

loOO  0.01 

The  values  of  AT  are  more  than  adequate  to  maintain  the 
accuracy  of  the  numerical  scheme.  The  minimum  number  of 
increments  for  a  natural  period  of  first  mode  lateral  vibra- 
tion of  the  column  in  the  numerical  process  Is  62.  The  values  of 
AT  were  selected  to  display  the  complete  history  of  the 
dynamic  response  to  the  pulse  loading  for  the  graphs  shown 
in  App.  III  using  the  maximum  900  points  permitted  in  Sub- 
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routine  GRAPH, 

The  values  of  maximum  stress  tabulated  on  the  following 
pages  are  maximum  compressive  stresses  if  there  is  no  sign 
before  the  number,  or  maximum  tensile  stresses  if  there  is 
a  negative  sign  before  the  number.   In  either  case  they  are 
the  maximum  absolute  values  of  stress  that  occurred  in  Spring 
1  or  Spring  2  of  the  column  model. 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c<-  2 


l/t>- 


1Q0 


e/^>«  0,05 


e 

Maximum 
Force 

<5n/*«  ) 

Residual 
Deflection 

(y»-/<?) 

Maximum 
Deflection 

(ym/9) 

Stress  1 
Maximum 

(s»max) 

Stress  2 
Maximum 

*  ^.rnaLX' 

0,01 

0.953 

2.12x10"* 

.  — '- 

0.00282 

0.00074 

0.05 

1.001 

2.24x10"* 

...  , 

0.00283 

0.00078 

0.10 

1.121 

2.50x10"' 

— — 

0.00287 

0.00100 

0.20 

1,319 

2.77x10"* 

-«.- 

0.00289 

-0.00129 

0.50 

1.658 

2.4.0x10'* 

— —  ■. 

0.00286 

-0.00233 

1.00 

1.935 

1.69x10"* 

— 

0.00278 

0.00213 

<=<-    2.5 


l/t>- 


100 


e/t>,    0.05 


? 

Maximum 
Force 

(?*/**  ) 

Residual 
Deflection 

(yr/<?) 

Maximum 

Deflection 

(y*n/£) 

Stress  1 
Maximum 

(stma.x) 

Stress  2 

■Maximum 

(sZmax) 

0.01 

0.953 

5.89x10'* 

•  — 

0.00312 

0.00084 

0.05 

1.029 

6.45x10"* 

-— 

0.00315 

-0.00099 

0.10 

1.189 

7.66x10"* 

a»«a 

0.00320 

-0.00132 

0.20 

1.470 

9.41X10"'* 

0.00327 

-0.00188 

0.50 

1.956 

8.88x10"* 

—  • 

0.00325 

r0.00275 

1.00 

2.331 

9.00x10"* 

0.00329 

-0.00263 

' 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


<*-  3 


t/<j-  100 


e/^„  0.05 


§ 

Maximum 
Force 

Residual 
Deflection 
(yr/9)  . 

Maximum 
Deflection 
(ym/9) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

0.01 

0.954 

1.14x10"' 

... 

0.00334 

-0.00110 

0.05 

1.069 

1.28x10"' 

... 

0.00337 

-0.00130 

0.10 

1.252 

1.54x10-' 

... 

0.00344 

-0.00167 

0.20 

1.602 

1.93x10"' 

... 

0.00352 

-0.00236 

0.50 

2.237 

2.20x10"' 

... 

0.00356 

-0.00310 

1.00 

2.633 

2.05x10*' 

...  : 

0.00360 

-0.00313 

<=<-  3.5 


ivba 


100 


e/^»  0.05 


? 

Maximum 
Force 

(P^/P*  ) 

Residual 
Deflection 

(yr/<?) 

Maximum 
Deflection 

(y*/?) 

Stress  1 
Maximum 

(stma.x) 

Stress  2 

-  Maximum 

(simax) 

0.01 

0.955 

1.81x10"' 

0.00349 

-0.00131 

0.05 

1.110 

2.05x10"' 

... 

0.00354 

-0.00155 

0.10 

1.314 

2.47x10"' 

... 

0.00360 

-0.00195 

0.20 

1.724 

3.25x10"' 

... 

0.00370 

-0.00268 

0.50 

2.447 

3.83x10"' 

MM 

0.00377 

-0.00340 

1.00 

2.6*51 

2.57x10"' 

... 

0.00380 

-0.00352 

• 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


• 


c<-  4 


Mj- 


100 


e/^«  0.05 


$ 

Maximum 
Force 

Residual 
Deflection 

(yr/9) 

Maximum 
Deflection 

(ym/9) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

'  <Lmax' 

0.01 

0.956 

2.54x10" 

»— ..  . 

0.003 61 

-0.00149 

0.05 

1.147 

2.89x10*' 

-— 

0.00366 

-0.00175 

0.10 

1.332 

3.49xlOH 

. 

0.00373 

-0.00219 

0.20 

1.636 

4.56x10-' 

— - 

0.003^3 

-0.00281 

0.50 

2.617 

5.57x10-' 

— - 

0.00392 

-0.00367 

1.00 

3.011 

2.93x10"' 

0.00394 

-0.00381 

ex.-  1 


i/t> 


70 


e/o*   0.05 


i 

Maximum 
Force 

(Pm/*.) 

Residual 
Deflection 

*  (yr/cO) 

Maximum 
Deflection 

(y»n/^) 

Stress  1 
Maximum 

(sima.x) 

Stress  2 
.Maximum 

(S^max) 

0.01 

0.841 

3.17x10"* 

0.00315 

0.00149 

0.05 

0.840 

3.19X10"2- 



0.00315 

0.00150 

0.10 

0.847 

3.13xlO~z 

--- 

0.00315 

0.00153 

0.20 

0.668 

4.13x10"* 



0.00324 

0.00163 

0.50 

0.934 

.2.83x10'* 



0.00312 

0.00190 

1.00 

0.962 

6.24xl0~3 

0.00272 

0.00224 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


<=<«  1.5 


1/^-  70 


e/^.  0.05 


$ 

Maximum 

Force 

(P^/P6  ) 

Residual 
Deflection 

Maximum 
Deflection 

(ym/^) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

( Szmax) 

0.01 

0.641 

2.5.0x10"' 



0.00361 

0.00149 

0.05 

0.654 

2.09x10"' 

... 

0.00361 

0.00153 

0.10 

0.69^ 

2.10x10"' 

0.00361 

0.00164 

0.20 

0.967 

2.10x10"' 

— 

0.00361 

0.00164 

0.50 

1.206 

2.04x10"' 

... 

0.00361 

-0.00260 

1.00 

1*337 

1. 03X10"1 



0.00364 

0.00314 

o< 


U<£ 


70 


eA»  0.05 


> 

Maximum 
Force 

(PW/PE  > 

Residual 

Deflection 

(yr/5) 

Maximum 
Deflection 

Stress  1 
Maximum 

(stma.x) 

Stress  2 

•Maximum 

( s2.max) 

0.01 

0.641 

3.79x10-' 

- — 

0.00404 

0.00149 

0.05 

0.676 

3.66x10-' 

wmmmwm^ 

0.00405 

0.00156 

0.10 

0.954 

3.96x10-' 

•»«•«» 

0.00406 

0.00176 

0.20 

1.060 

4.05x10"' 

«»a»«» 

0.00407 

-0.00226 

0.50 

1.377 

3.43x10"' 

... 

0.00403 

-0.00366 

1.00 

1.516 

2.22x10"' 



0.00402 

0.00364 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c<-  2.5 


!/<£-  70 


e/(3„  0.05 


■ 

Maximum 
Force 

(P~/PK  ) 

Residual  . 
Deflection 

Maximum 
Deflection 
(ytn/t^) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

0.01 

0.S41 

5.36x10"' 

0.00418 

0.00149 

0.05 

0.397 

i>.52xlO"'  . 

— 

0.00419 

0.00161 

0.10 

0.995 

5.77x10"' 

-.— 

0.00421 

-0.00202 

0.20 

1.167 

5.99xlOH 

1.67 

0.00423 

-0„ 00308 

0.50 

1.46*9 

5.03xlOH 

... 

0.00422 

-0.00413 

1.00 

1.619 

3.06xKT' 

— 

0.00422 

0.00393 

.' 


c<-  1 


lyfc- 


50 


e/^,  0.05 


s 

Maximum 
Force 

<P«/*.> 

Residual 
Deflection 

(yr/to). 

Maximum 
Deflection 

(y»n/0 

Stress  1 
Maximum 

Stress  2 

■Maximum 

(s2.max) 

0.01 

0.585 

2.67x10"' 



0.00417 

0.00234 

0.05 

0.592 

2.69x10'' 

... 

0.00418 

0.00240 

0.10 

0.614 

2.70xl0H 

... 

0.00419 

0.00251 

0.20 

0.655 

2.63x10"' 

— 

0.00417 

0.00275 

0.50 

0.731  • 

2.52x10"' 

™ 

0.00419 

0.00329 

1.00 

0.774 

8.60xlO"1 

0.00396 

0.00366 

- 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c<-    1.5 


1/<J- 


50 


e/^,    0.05 


" 

Maximum 

Force 

(P^/P6  ) 

Residual 
Deflection 

(y*-/<?) 

Maximum 
Deflection 

(ym/9) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

0,01 

0.535 

5.17x10"' 

.__ 

0.00446 

0.00234 

0.05 

0.612 

5. 16x10-' 

.  — 

0.00446 

0.00251 

0.10 

0.654 

'5.13x10-' 

'."*■" 

0.00446 

0.00275 

0.20 

0.708 

5.05x10"' 

—  . 

0.00446 

0.00307 

0.50 

0.822 

4.03x10"' 

... 

0.00445 

0.00383 

1.00 

0.864 

1.60X10-1 



0.00435 

0.00413 

c< 


I/to- 


50 


e/o»     0.05 


? 

Maximum 
Force 

Residual 
Deflection 

(yr/co) 

Maximum 
Deflection 

(y*>/£) 

Stress  1 
Maximum 

Stress  2 

Maximum 

(sZmax) 

0.01 

0.585 

7.27x10"' 

... 

0.00461 

0.00235 

0.05 

0.626 

7.34x10" 



0.00462 

0.00257 

0.10 

0.679 

7.37x10"' 



0.00462 

0.00287 

0.20 

0.756 

7.18x10"' 

0.00462 

0.00338 

0.50 

0.873 

4.93x10"' 



0.00461 

-0.00454 

1.00 

0.913 

2.07x10"' 



0.00456 

0.00439 

V-8 
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APPENDIX  V   . 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


<=<-  2.5 


1/<J-  50 


e/^„  0.05 


§ 

Maximum 
Force 

(Pm/P8  ) 

Residual 
Deflection 

(yr/9) 

Maximum 
> Deflection 

(y*>/c>) 

Stress  1 
Maximum 

(simax' 

Stress  2 
Maximum 

0.01 

0.586 

9.14x10"' 

— 

0.00472 

0.00236 

0.05 

0.645 

9.30x10"' 

— 

0.00473 

0c 00269 

0.10 

0.700 

9.40x10"' 

— 

0.00473 

0.00300 

0.20 

0.794 

8.98x10"' 

0.00474 

-0.00413 

0.50 

0.909 

6.17x10-' 

0.00475 

-0.00475 

1.00 

0.948 

2.41x10"' 

0.00471 

0.00456 

o<-  2 


1/^. 


100 


e/^.  0.01 


$ 

Maximum 
Force 

(P*/P*  ) 

Residual 
Deflection 

(yr/<?) 

Maximum 
Deflection 

(y»n/£) 

Stress  1 
Maximum 

(simax) 

Stress  2 

•  Maximum 

(s2.max) 

0.01 

1.003 

2. 56x10"* 



0.00287 

0-.  00100 

0.05 

1.192 

3.03x10"* 

— . 

0.00292 

0.00118 

0.10 

1.376 

3.29xlO"1 

1.64 

0.00295 

0.00149 

0.20 

1.616 

3.30x10"* 

1.81 

0.00295 

-0.00180 

0.50 

1.911 

3.07xlO-1 

1.93 

0.00294 

-0.00231 

1.00 

1.988  . 

2. 95x10" 3 

0.69 

0,00250 

0.00235 

V-9 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


<=*-  2.5 


L/e>-  100 


e/^«  0.01 


e 

Maximum 
Force 

Residual 
Deflection 

(y*-/<?). 

Maximum 
Deflection 

(ym/9) 

Stress  1 
Maximum 

(simax' 

Stress  2 
Maximum ' 

'  ^.max' 

0.01 

1.017 

6.93xl0"2 

— 

0.00317 

0.00100 

0.0$ 

1.256 

8.89X10"2" 

mm  mm,  ^m 

0.00325 

-0.00139 

0.10 

1.486 

1.04x10-' 

0.00331 

-0.00185 

0.20 

1.813 

1.17x10-' 

... 

0.00334 

-0.00250 

0.50 

2.264 

1.11x10-* 

... 

0.00335 

-0.00286 

1.00 

2.423 

2.29x10"* 

0.00306 

0.00287 

c<«  3 


V<£ 


100 


e/^»  0.01 


i 

Maximum 
Force 

Residual 
Deflection 

(yr/^) 

Maximum 
Deflection 

(y*n/£) 

Stress  1 
Maximum 

Stress  2 

■Maximum 

(szmax) 

0.01 

1.032 

1.30x10"' 

... 

0.00338 

-0.00115 

0.05 

1.320 

1.72x10"' 

... 

0.00348 

-0.00173 

0.10 

1.578 

2.07x10"' 

... 

0.00354 

-0.00224 

0.20 

1.971 

2.47x10"' 

... 

0.00360 

-0.00274 

0.50 

2.535 

2.11xl0H 

.... 

0.00359 

-0.00326 

1.00 

2.744 

7.27xl0"»- 

0.00346 

0.00327 

V-10 


..  .' 


APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


<*-  3.5 


I/t> 


100 


e/^o.  0.01 


§ 

Maximum 
Force 

(P~/Ps  ) 

Residual 
Deflection 

Maximum 
Deflection 
(y.n/9) 

Stress  1 
Maximum 

(sirr»ax) 

Stress  2 
Maximum 

0.01 

1.046 

2.00x10"' 

0.00353 

-0.00135 

0.05 

1.360 

2.69x10"' 

— 

0.00363 

-0.00199 

0.10 

1.666 

3.26xlO'j 

:  — — 

0.00370 

-0.00255 

0.20 

2.106 

3.66x10"' 

0.00377 

-0.00292 

0.50 

2.739 

3.00x10"' 

— 

0.00376 

-0.00356 

1.00 

2.961 

1.35x10-' 

-.— 

0.00372 

0.00354 

c< 


i/q* 


100 


e/^= 


0.01 


? 

Maximum 

Force 

(P^/P*  ) 

Residual 
Deflection 

Maximum 

Deflection 

(y»n/£) 

Stress  1 
Maximum 

Stress  2 

■Maximum 

( Sz.max) 

0.01 

1.062 

2.74x10"' 

--- 

0.00364 

-0.00151 

0.05 

1.433 

3.70x10-' 

... 

0.00375 

-0.00220 

0.10 

1.751 

4.52x10-' 

2.67 

0.003 S3 

-0.00273 

0.20 

2.230 

5.49x10"' 

2.97 

0.00390 

-0.00306 

0.50 

2.694 

3.94x10-' 

3.54 

0.00366 

-0.00363 

1.00 

3.114 

1.93X10-*1 

2.70 

0.00390 

0.00374 

V-ll 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c*.  1 


l/tj-  70 


e/^«  0,01 


§ 

Maximum 
Force 

<P„/P6  ) . 

Residual 
Deflection 

Maximum 
Deflection 

(ym/9) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

0.01 

0.916 

2.54xlO-2- 



0.00309 

0.00197 

0.05 

0.927 

2. 89x10" * 

«»^«» 

0.00313 

0.00201 

0.10 

0.947 

3.63x10."* 

0.49 

0.00320 

0.00206 

0.20 

0.973 

2.28x10"* 

0.45 

0.00306 

0.00217 

0.50 

0.991 

3.94xl0"3 

0.20 

0.00264 

0.00232 

1.00 

0.994 

Iil3xl0"3 

0.06 

0.00252 

0.00242 

<=<-  1.5 


L/^-  70 


»/^»   0.01 


. 


$ 

Maximum 
Force 

(JU/p*  ) 

Residual 
Deflection 

(yr/^) 

Maximum 
Deflection 

Stress  1 
Maximum 

Stress  2 
■Maximum 

( s2.max) 

0.01 

0.922 

2.24x10"' 



0.00383 

0.00199 

0.05 

0.994 

2.24x10-' 



0.00383 

0.00212 

0.10 

1.072 

2.22x10-' 



0.00383 

0.00230 

0.20 

1.181 

2.17x10"' 



0.00383 

0.00261 

0.50 

1.333 

1.60x10'' 

—  . 

0.00377 

0.00315 

1.00 

1.364 

2. 11x10" l 



0.00345 

0.00334 

V-12 
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APPENDIX  V 
Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


C<-  2 


l/t>-  70 


e/^a  0.01 


e 

Maximum 
Force 

Residual 
Deflection 

Maximum 
Deflection 

(ym/9) 

Stress  1 
Maximum 

(stmax) 

Stress  2 
Maximum 

0.01 

0.930 

4.01x10"' 

... 

0.00406 

0.00200 

0.05 

1.041 

4.17x10"' 



0.00406 

0.00222 

0.10 

1.151 

4.17x10"' 

mm  «  •» 

0.00408 

0.00250 

0.20 

1.296 

4.02x10"' 

0.00408 

-0.00309 

0.50 

1.497 

3.04x10-' 

. 

0.00408 

0.00359 

1.00 

1.542 

5.20xlO"1 



0.00388 

0.003 80 

c<«  2.5 


Iv^« 


7Q 


e/k»  0.01 


.■?: 

Maximum 
Force 

Residual 
Deflection 

(yr/c^) 

Maximum 
Deflection 

(y*n/?) 

Stress  1 
Maximum 

(Sim;ix) 

Stress  2 
-  Maximum 

( s2.max) 

0.01 

0.939 

5. 62x10"' 

0.00420 

0.00201 

0.05 

1.079 

6.02x10"' 

• 

0.00423 

0.00231 

0.10 

1.210 

6.17x10*-' 

1.85 

0.00424 

-0.00292 

0.20 

1.382 

6.03x10"' 

1.95 

0.00425 

-0.00364 

0.50 

1.594 

3.95x10"' 

1.99 

0.00425 

-0.00396 

1.00 

1.642 

7.93x10*-* 

0.72 

0.00413 

0.00405 

v-13 


Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


l/% 


50 


V^; 


0.01 


'■ 

Maximum 
Force 

(?~A  ) 

Residual 
Deflection 

Maximum 
Deflection 

Stress  1 
Maximum 

(sinnax) 

Stress  2 
Maximum 

*  ^.rnaix' 

0.01 

0.633 

2. 71x1 0"' 



0.00419   0.00284 

0.05   0.658 

2.74x10"' 

-_- 

0.00420 

0.00297 

0.10 

0.689 

2.75x10-' 

0.67 

0.00421 

0.00313 

0.20 

0.730 

2.65x10"' 

0.75 

0.00421 

0.00337 

0.50 

0.775 

I.l5xl0-1 

0.48 

0.00403 

0.00367 

1.00   0.782 

1.78x10-* 

0.10 

0.00384 

0.00378 

o< 


1.5 


1/^- 


50 


e/p 


0.01 


$ 

Maximum  j  Residual 
Force   Deflection 

Maximum 
Deflection 

(ym/^) 

Stress  1 

Maximum 

Stress  2 

-  Maximum 

(s2.max) 

0.01 

O.636    5.30x10" 



0.00447 

0.00286 

0.05 

0.687    5.22x10-' 



0.00447 

0.00311 

0.10   0.735 

5.02x10"' 

1.12 

0.00447 

0.00338 

0,20 

0.791 

4.57x10-' 

1.10 

0.00447 

0.00369 

0.50 

0.859 

2.62x10-' 

0.90 

0.00443 

0.00411 

1.C0 

0.870 

3.66x10-* 

0.16 

0.00426 

0.00422 

.  V-14 


•APPENDIX  V 
Maximum  Force,   Residual  Deflection  and  Maximum  Stress  Data 


OC-    1 


l/(j- 


50 


e/^„     0 


.01 


§ 

Maximum 
Force 

Residual 

Deflection 

(y*-./<?) 

Maximum 
f Deflection 

(ym/9) 

Stress  1 
Maximum 

(simax) 

Stress  2 
Maximum 

0.01 

0.633 

2.71x10"' 



0.00419 

0.00234 

0.05 

0.653 

2.74x10"' 

0.00420 

0.00297 

0.10 

0.639 

2.75X10"1 

0.67 

0.00421 

0.00313 

0.20 

0.730 

2.65x10"' 

0.75 

0.00421 

0.00337 

0.50 

0.775 

1.15xl0-' 

0.43 

.0.00403 

0.00367 

1.00 

0.732 

li 73x10" * 

0.10 

0.00334 

0.00373 

. 


o<-     1.5 


I^> 


50 


e/^     0 


.01 


i 

Maximum 
Force 

Residual 

Deflection 

(yr/5) 

Maximum 
Deflection 

(y»n/£) 

Stress  1 

Maximum 

Stress  2 

-  Maximum 

(s2.max) 

0.01 

0.636 

5.30x10"' 



0.00447 

0.00236 

0.05 

0.637 

5.22x10"' 

0.00447 

0.00311 

0.10 

0.735 

5.02x10"' 

1.12 

0.00447 

0.00333 

0.20 

0.791 

4.57x10"' 

1.10 

0.00447 

0.00369 

0.50 

0.359 

2.62x10-' 

0.90 

•  0.00443 

0.00411 

1.00 

0.370 

3.66x10"* 

0.16 

0.00426 

0.00422 

.  V-14 


Maximum  Force, 


APPENDIX  V 
Residual  Deflection  and  Maximum  Stress  Data 


OC-  2 


L/^-  50 


e/^«  0,01 


§ 

Maximum 
Force 

Residual  1 
Deflection 

(yr/<?) 

Maximum 
Deflection 

(ym/c^) 

Stress  1 
Maximum 

(simax) 

Stress  2 

Maximum 

0.01 

0.641 

7.44xl0H 

— - 

0.00463 

0.00289 

0.05 

0.709 

7.44x10-' 



0.00463 

0.00323 

0.10 

0.766 

7.17x10"' 



0.00464 

0.00354 

0.20 

0.632 

6.27x10-' 

... 

0.00463 

0.00391 

0.50 

0.907 

3.44x10"' 



0.00463 

0.00435 

1.00 

0.919 

4. 44x10" t 



0.00449 

0.00446 

<=<-  2.5 


l/^»  50 


e/o=.  0.01 


e 

Maximum 
Force 

Residual 
Deflection 

(yr/<?) 

Maximum 
Deflection 

(y*n/?) 

Stress  1 
Maximum 

(stm;a.x) 

Stress  2 

. Maximum 

( s2.max) 

0.01 

0.646 

9. 35x10"' 



0.00473 

0.00291 

0.05 

0.729 

9.45x10-' 

0.00475 

0.00334 

0.10 

0.7^9  . 

9.19x10"' 

1.65 

0.00476 

0.00367 

0.20 

0.863 

7.94x10"' 

1.60 

0.00476 

-0.00427 

0.50 

0.941 

4.03x10"' 

1.26 

0.00476 

0.00452 

1.00 

0.952 

5.25x10-*- 

0.22 

0.00465 

0.00462  | 
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APPENDIX  VI 
Evaluation  of  the  Approximations  and  of  the  Accuracy  of 
the  Numerical  Scheme 

From  a  consideration  of  the  maximum  lateral  deflections 
encountered  in  the  solutions  the  approximations  made  in  de- 
veloping the  equations  for  the  column  model  may  be  evaluated. 
The  maximum  value  of-©-  for  all  solutions  was  less  than  0.08 
radian.   Based  on  a  maximum  -©■  =  0o08  radian  the  errors  are 
as  follows: 

(a)  0„01  per  cent  error  in  (L/2-e  tan0  )  as  a  result 
of  neglecting  e  tan-O- 

(b)  0.2  per  cent  error  in  |-^-sin-©,-e(l-cos-©')J  by  con- 
sidering sln-©-=-©-  and  neglecting  e(l-cos-O) 

(c)  0.25  per  cent  error  in  tanOand  0.11  per  cent 
error  In  sin-O1  as  a  result  of  considering 
tan  -©■  =sin-0-=*©- 

(d)  less  than  0.06  per  cent  error  in  (l-cos-0-)  as  a 
result  of  considering  l-cos^=-£-0* 

From  this  the  approximations  which  have  been  made  in  deriv- 
ing the  equations  for  the  column  model  are  reasonable  and 
produce  insignificant  error  in  the  results. 

As  an  additional  check  for  cumulative  error  a  compar- 
ison of  the  ram  work  done  on  the  column  model  with  the  net 
energy  input  to  each  of  the  springs  (a  combination  of  strain 
energy  of  plastic  deformation  and  the  existing  elastic  strain 
energy)  is  made  in  all  computer  solutions  for  the  results 
shown  in  Figs.  10  thru  14„   These  agree  within  1  per  cent 
thus  verifying  the  compatibility  of  the  equations  and  check- 


VI-I 


ins  tne  accuracy  of  the  numerical  scheme  for  the  time  incre- 
ments chosen.   Work-energy  considerations  are  not  included 
in  the  computer  program  for  obtaining  graphical  results  but 
the  dimensionless  time  increment  AT  used  for  the  numerical 
integration  in  both  programs  for  a  given  value  of  S  was  the 
same. 
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